Cusps and shocks in the renormahzed potential of glassy random manifolds: How 
Functional Renormalization Group and Replica Symmetry Breaking fit together 
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We compute the Functional Renormalization Group (FRG) disorder-correlator function R{v) for 
d-dimensional elastic manifolds pinned by a random potential in the limit of infinite embedding 
space dimension A''. It measures the equilibrium response of the manifold in a quadratic potential 
well as the center of the well is varied from to ii. We find two distinct scaling regimes: (i) a "single 
shock" regime, v"^ ~ L""* where L"^ is the system volume and (ii) a "thermodynamic" regime, 
V r^ N . In regime (i) all the equivalent replica symmetry breaking (RSB) saddle points within the 
Gaussian variational approximation contribute, while in regime (ii) the effect of RSB enters only 
through a single anomaly. When the RSB is continuous (e.g., for short-range disorder, in dimension 
2 < d < 4), we prove that regime (ii) yields the large- A*' FRG function obtained previously. In 
that case, the disorder correlator exhibits a cusp in both regimes, though with different amplitudes 
and of different physical origin. When the RSB solution is 1-step and non-marginal (e.g., d < 2 
for SR disorder), the correlator R{v) in regime (ii) is considerably reduced, and exhibits no cusp. 
Solutions of the FRG fiow corresponding to non-equilibrium states are discussed as well. In all cases 
the regime (i) exhibits a cusp non-analyticity at T = 0, whose form and thermal rounding at finite 
T is obtained exactly and interpreted in terms of shocks. The results are compared with previous 
work, and consequences for manifolds at finite A*', as well as extensions to spin glasses and related 
models are discussed. 
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I. INTRODUCTION 

A major difficulty in devising analytical methods to 
handle glassy systems, such as systems with quenched 
disorder, is to describe accurately the many metastable 
states which play a role both in the statics (equilibrium) 
and the dynamics, as well as the barriers separating 
these states. Two main general methods have been 
developed. The first is a mean-field theory, based on the 
Gaussian variational method (GVM), which, in the stat- 
ics, captures the many states by a multitude of saddle 
points exhibiting spontaneous replica symmetry break- 
ing (RSB)i. Being highly versatile, the GVM has been 
applied to numerous problems, notably in spin glasses, 
and has later been extended to the dynamics^.!.?.. The 
second method is the functional renormalization group 
(FRG) which has been applied successfully to disordered 

elastic systems and random field spin models, both in the 
staticsi>^>^>^>^>^ili>^>H^>i^as^^>i^>20^^^^^^^, 

as for the driven dynamics2^i2£i2iiHi^i2ii2^i2&. Although 
its range of applications is at present smaller, it is a 
powerful and promising method which allows to compute 
fluctuations not captured by mean-field theory. The 
(relevant) coupling constant in FRG is the disorder 
correlator R{u). In contrast to standard field theories 
it is not a number but a function of the field u, defined 
for the microscopic model in Eq. ([5]). Since the two 
methods GVM and FRG are rather different in spirit, 
and historically have developed along separate tracks, 
it is important to compare them whenever possible. A 
further goal is to understand whether and how they may 



be extended to a broader range of models. 

In this paper we focus on disordered elastic systems, 
where both methods have been applied. As we discuss, 
some of the conclusions and ideas may extend to other 
models with quenched disorder. Besides being of direct 
interest for experiments, including vortex lattices, mag- 
netic systems, density waves^ i'^^i'^^'^^i^°'^^i'^^ , models of 
manifolds in a random potential provide the simplest 
example of a glass phase where numerous metastable 
states occur beyond the so-called Larkin scale. As for 
random-field systems, the so-called dimensional reduc- 
tion phenomenon occurs, which renders conventional 
zero-temperature perturbation theory trivial^i^i^, indi- 
cating the failure of the latter to capture the complexity 
of the energy landscape. A big success of both methods, 
GVM and FRG, has been to circumvent this problem. 
However, they achieve this in seemingly rather different 
ways: The GVM approximates the Gibbs measure by a 
hierarchical superposition of Gaussians, encoded in the 
Parisi Ansatz^''*^. In the FRG, the existence of a cusp 
(in the coupling function R" {u)) beyond the Larkin scale 
is related to the existence of many metastable states. 
This was nicely illustrated in an early paper by Balents, 
Bouchaud and Mezard^. The idea that coarse graining 
leads to shocks in the force landscape, and the similarity 
of the pinning problem to the Burgers equation, was in- 
troduced using a toy RG model with two degrees of free- 
dom. However, the function R(u) defined in that work 
does not coincide with the one usually studied in field 
theory, making a precise comparison difficult. 

In this article, we want to make this comparison quan- 



titative. Let us denote u{x), x G R"^, the A''-component 
displacement (or "height" ) field which parameterizes the 
position of a manifold of internal dimension d in the iV- 
dimensional embedding space. The GVM was applied 
to the problem by Mezard and Parisi (MP)^^, introduc- 
ing replicated fields Ua{x). For long-range disorder, or 
for short-range disorder and internal dimensions above 
d = 2, they found a solution with continuous replica sym- 
metry breaking (RSB) with a roughness exponent for the 
manifold u{x) ~ x^ given by the Flory estimate. (The 
FRG allows one to go beyond this result and to com- 
pute deviations from Flory for A'^ < oo.) For short-range 
disorder, with d < 2, they found a 1-step RSB solution, 
analogous to the one in infinite range p-spin modela^*^. 

In the infinite- A^ limit, the GVM becomes formally ex- 
act and hence can be compared with the FRG. Two of 
usiiiH obtained a self-consistent equation for the cou- 
pling function of the FRG, R{u). It was defined from the 
effective action of the replicated field theory (for a uni- 
form mode u{x) — u), a standard field theory definition, 
and computed in the large- A'^ limit, performing the usual 
rescaling R{u) = NB{u'^/N). Although the resulting 
self-consistent equation for B{x) is formally valid only 
below the Larkin scale, the corresponding FRG equa- 
tion could be continued "naturally" to scales beyond the 
Larkin scale (with a cusp present at u = 0), extend- 
ing the flow all the way to the RG fixed point. For 
2 < d < 4, where the RSB is continuous, this FRG flow 
recoverediiii^ the MP result for small overlap, i.e., it 
yielded only, even though exactly, the (non-trivial) con- 
tribution of the most distant states to the correlation 
function. However, quite surprisingly, varying the IR cut- 
off m in the confining potential well finally allows for the 
reconstruction of the complete self-energy function ob- 
tained in the GVM, without ever referring to ultrametric 
matriceSiiiii^ . 

Despite this quantitative progress in understanding the 
connection between FRG and GVM, several questions re- 
mained. It is natural that the FRG recovers the contribu- 
tion of distant states, since it introduces an external field 
which explicitly breaks replica symmetry and hence splits 
all the replicas: Ua — Ub ^ 0. However, one would hope 
FRG to describe all states, not just the ones with smallest 
overlap. Further, no thermal rounding of the cusp was 
found which is physically surprising in view of results re- 
lating finite-temperature droplets and FRG^i^i^. Fi- 
nally, one would like to describe better the situation 
where the GVM yields a non-marginal 1-step RSB so- 
lution, namely the case d < 2 with short-range disor- 
der, which includes in particular the KPZ problem with 

d = x50-51.52.53.54 'j 

The first aim of this paper is to compute the FRG 
functions from first principles in the large- A^ limit. We 
take advantage of the recently obtained direct relation 
between the field theoretic definition of the FRG function 
R{u) and directly observable quantities'**''*^. This has 
allowed for a numerical determination of R{u) or the force 



correlator 



A(u) = -R"{u) 



(1) 



for A^ = 1 interfaces at T = in dimensions d = 0, 1, 2, 3 
in Ref. 55, as well as for the dcpinning problem in 
Refs. 56.57 . It was found that the numerics compares 
remarkably well with the e — A — d expansion, to 1-loop, 
and even better to 2-loop order. The idea, also imple- 
mented here, is to subject the manifold, in addition to 
the random potential, to an external quadratic potential 



well ^ ^ d X (u[x) — w)^, centered at v. The "mass" m 
acts as an infrared cutoff limiting the interface fluctua- 
tions. By measuring the free energy V{v) of the system 
as a function of v one obtains a random landscape whose 
second cumulant is the function R(v). More generally, if 
V — > v{x) the whole second cumulant functional R[v] is 
retrieved. 

Here we compute this functional exactly. We find two 
distinct scaling regimes with non-trivial infinite- A^ limits, 
the first one where w^ ^ N^ L~'^, the second for w^ ^ A^. 
The reason for this peculiar property is that at t; = 
there is spontaneous RSB, which implies the contribution 
of many saddle points which are equivalent under replica 
permutations. If the "applied field" v remains "small" 
(first regime) the spontaneous RSB saddle point is not 
significantly modified and all saddle points contribute to 
a given observable (though, now, not all of them equiv- 
alently) . This can be handled by a method introduced 
in Ref. |43 , and will be applied here with some improve- 
ments. In the second regime, which corresponds to the 
more conventional scaling at large A^, the applied field is 
stronger and the saddle point is modified. We explicitly 
compute the w-dependence, and show how, for large v, 
a non-trivial FRG function emerges. Specifically, in the 
case of a uniform field v{x) = u we obtain: 



R{v) - i?(0) 



L-'iriv^L'^) for v^ 



L- 



(*) 



Nr{v^/N) for v^ - A^, (m) 



(2) 



with two different scaling functions which we compute. 
We check that the two regimes match, i.e., f{z) ~ Az at 
large z and r{z) ^ Az at small z. In the case where the 
RSB is continuous (e.g., for 2 < d < 4) we prove that 
regime {ii) yields the large- A^ FRG function obtained in 
our previous stud y^^i^^ . Hence the natural continuation 
(using the FRG flow-equation) performed there is cor- 
rect, and one of the main results of the present paper is 
to show this rigorously. Remarkably, RSB enters in this 
regime only through a single number, an anomaly. In 
the case of continuous RSB (including the marginal 1- 
step solution in d = 2) we find that both regimes (i) and 
[ii] in ([2]) exhibit a cusp in R"{v) at the origin, though of 
different nature. Specifically, one finds f{z) = Bz+Cz'^/'^ 
at small z (at T = 0) and r{z) ^ Az + Cz^l"^ at small z 
(for all T in the glass phase). These cusps are caused by 
jump discontinuities in V'iv), called shocks, as the center 
of the well v is moved. In regime (j) the scaling function 
describes a fiuctuation of the energy, as measured by the 



connected correlator 



[V{v) - l/(0)]2 = 2L''[R{v) - R{0)] 



(3) 



of order unity as v is moved by a very small amount ~ 
j^-d/2^ Hence we call this regime "single shock regime". 
By contrast, in regime [ii) the variance of the energy fluc- 
tuations scales with volume and N . Hence we call this 
regime "thermodynamic" , since the scaling function r{z) 
encodes the average properties of many shocks. The two 
regimes also exhibit different properties with respect to 
temperature. In regime {it), as was found i n^^i^'^ , the 
temperature dependence is weak and the cusp survives 
even for T > 0. Regime (i) exhibits a thermal rounding of 
the cusp, whose form is obtained in an exact closed form. 
We compare its form with the predictions for the ther- 
mal boundary layer of the FRG obtained previously^i^ 
in cases where it can be obtained using droplet arguments 
(i.e., at finite N). For sufficiently short-ranged disorder 
in d < 2, the equilibrium 1-step RSB solution yields a 
considerably reduced regime (m) with no non-analyticity. 
A non-analyticity in this regime would indicate the crit- 
icality of the system towards clustering of the states into 
an ultrametric superstructure (higher-step RSB). This 
kind of criticality is not present in typical samples in 
d < 2, and only arise in exponentially rare disorder re- 
alizations, as we confirm by studying the configurational 
entropy of non-equilibrium states. Possible consequences 
of our results for manifolds at finite N are discussed, as 
well as extensions to spin glasses and related models. 

The outline of the paper is as follows. In Section |TT] 
we define the model and the observable to be studied. 
In Section IIIII we analyze regime (i) for both continu- 
ous and 1-step RSB, discussing various subtleties of the 
phase diagram in the 1-step case. We study in particu- 
lar the non-analytic cusp arising in the force correlator 
at r = and its thermal rounding. Taking the limit of 
large w ^ 1 we establish the matching with the regime 
v^ ^ N, which is analyzed in detail in Section IIVI We 
rigorously show how to compute the FRG function R{v) 
exactly by introducing two replica groups and derive from 
it the correct FRG flow equations in all cases, including 
the anomaly arising from RSB. The physical significance 
of the presence or absence of a cusp in this regime is 
discussed. The results are summarized in Section IVl and 
possible applications are discussed. 



II. MODEL, OBSERVABLES AND PREVIOUS 
RESULTS 



A. Model 

We consider an elastic manifold parameterized by a 
A''-component displacement field u{x), also denoted u^, 
where x belongs to the internal d-dimensional space. 
The manifold is exposed to a random potential, V{x, u), 
which lives in a (D = d + A^)-dimensional space. Indices 
of the field u]., i — I, . . . , N are shown only when strictly 



needed, and we use the notation u ■ v = X]i=i ■"'^'- We 
study the (classical) equilibrium statistical mechanics de- 
fined by the canonical partition sum Zy = Tr e^^^^' at 
temperature T, and denote thermal averages by {F[u\)v 
(or sometimes simply {F[u\)) in a given realization of 
the random potential. The model is defined by the total 
energy: 



Hv[u] 



g^ ^Uk ■ u-fc + 



V{x,Ux) , 



(4) 



^—ikx 



where Uk = J^Uxe''^^, J^ = J d''-x and Ux = Jk'^ke' 
Ik ^ / (27r)'' • "^^ ^-^ *^^ average center-of-mass position 
Ucm '■— Uk=o/L'^ we choose a non-zero value for g^^Q — 
cm?, which takes the role of a mass, c being the elastic 
constant. The mass provides a quadratic well for the 
manifold and thus serves as an IR cutoff to limit the 
displacement fluctuations. One is often interested in the 
scale invariant limit, r/i — > 0. A UV cutoff ~ A~^ in x 
space (e.g., due to a lattice) is implicit everywhere. For 
specific applications we consider— : 



1 



9k 



c{k^ + m?) 



(5) 



even though most results apply to more general forms of 
gk ■ The quenched disorder is chosen to possess statistical 
translational invariance, with second cumulant 



V{x, u)V{x', u') = S^'^^ {x ~ x')Ro{u - u') 



(6) 



This property entails a useful symmetry - see below - 
usually referred to as statistical tilt symmetry (STS). We 
always assume 0{N) symmetry of the disorder, choosing 
the bare correlator to be^: 



Ro{u) = NB{u^/N). 



(7) 



This scaling with N yields a non-trivial large- A^ limit. 
Among the variety of models parameterized by the func- 
tion B{z) one distinguishes short-range (SR) disorder, 
one often studied example being: 



(I) Biz) ^ Boe-'/"-" [SR], 



(8) 



and long-range (LR) disorder, often represented by the 
family of power law force correlator—: 



(II) B'{z) 



Bn 



1(1 



"fV 



.)7 



[LR]. 



(9) 



Here, the parametrization of model II is chosen^** such 
that the limit 7 — > cxd corresponds to model I, and that 
the limits 7^1 and 7^0 are meaningful. These 
two models possess a special scale-invariance property at 
infinite A^: As discussed below in Section HVl they arrive 
at their FRG fixed point after a finite renormalization 



time 



85 



It is well known— that the effect of disorder in model 
(m for any N becomes non-linear, and metastability ap- 
pears when the mass is decreased beyond some charac- 
teristic scale. As easily seen from dimensional analysis, 
the natural unit is the so-called Larkin scaled: 



1 

Mc = ^ := 






lA 



(10) 



where e = A — d. In finite dimensions d > 0, Lc = I/Mc 
has the loose interpretation of the smallest typical size of 
domains which may be trapped in different metastable 
states and thus exhibit glassiness at low temperature. 
The energy of such domains is naturally expressed in the 
unit of energy 



Er 



cfj-c '^ (rft^cy 



(11) 



while rf and Lc — 1/ fJ-c sue the natural scales for embed- 
ding space and internal space (i.e. inverse mass) lengths 
in the problem. We are free to choose units in which 
Ec = Lc = Tf ^ 1, ov equivalently, c = Bq = Vf = 1, 
which we adopt throughout the paper. For complete- 
ness, we give the dimensions of all observables used in 
the present paper in App. \^ which allows to restore the 
full dependence on these parameters. 

In order to study the model ([4]) one introduces repli- 
cated fields u°'{x), a — 1, . . . ,n. Using standard meth- 
ods, all disorder-averaged correlation functions of the 
u{x) field can be expressed as correlation functions of 
the replicated fields u°'{x) in the theory with partition 
sum Zy = Tre""^!"! and action S[u\: 



S[u]^ 



rp 2-^ 



9l'K 



1 

2r2 



J2 f RoK 



(12) 



B. Summary of previous studies 

1. Gaussian Variational Method (GVM) 

Before defining the observable computed in the present 
work let us briefly review the quantities studied in pre- 
vious publicationsiiii^i^, and the main results obtained 
there (details are skipped and can be found in these origi- 
nal publications). The model p^ was studied in Ref. 149 
using the Gaussian Variational Method (GVM), which 
becomes exact in the limit N = oo. The central observ- 
able calculated there is the two point correlation function 
of the replicated field: 

(u^'4^)s - Gabim,{2TTf6^''\k + k'). (13) 

Here we denote by {.)s averages over the replicated ac- 
tion (we later drop the subscript S when not strictly 
needed). The correlation P^ encodes^ the averages 

{ulvP,^,) = Gaa{k)5,j{2'KY5('^'>{k + k') (diagonal part) 



and {ului.Y = Y.bGab{k)5ij{2TTY5^'^\k + k') (connected 
thermal average) . The large distance behavior of the first 
one defines the roughness exponent C, of the manifold, i.e. 
((w^-02)~ |x-x'pC, equivalently Gaa(yfc) - fc-(''+20 
at small k. These hold at fixed scale in the limit to ^ 0, 
or at small but fixed m for scales smaller than the IR 
cutoff, k 3> 1/m. Another important exponent charac- 
terizes how the fluctuations of the ground state energy 
(or of the free energy) scale with system size, /S.E ~ L^ . 
Here, thanks to the statistical tilt symmetry, one has the 
relation 9 = d - 2 + 2^. 

Quite generally P^ takes the form: 



TGjik) = Sab9k ^ - <^ab{k) , 



(14) 



and within the GVM the self-energy is taken to be k- 
independent, aab{k) = Uab- It obeys a self-consistent 
equation which arises as a large- iV saddle-point equation, 
reading: 



^ab = -7fB' 



[Gaa{k)~Gab{k)] 



(15) 



There are two types of saddle points: Either they re- 
spect the replica symmetry of the action p2|) . Oa^b = cr, 
which happens in the high-temperature phase (where the 
roughness exponent assumes its thermal value, C = Cth = 
max(0, (2 — d)/2). Or, the saddle points spontaneously 
break the symmetry (RSB), aab -^ o'(u), where u G [0, 1] 
labels the distance of replicas in an ultramctric Parisi 
scheme describing the glassy, pinned phase. 



Let us summarize the results for N — oo and 



e 



9k = 



for which the dependence on the mass is worked 



out in Ref. 13. 

We start with the case where the glass phase is de- 
scribed by continuous RSB (also called infinite-step or 
full RSB). Within such a RSB scheme it is found that 
the roughness exponent equals its "Flory" value: 



C = Cf = 



(I) 

(4-d)/[2(l + 7)], (II) 



(16) 



and the continuous RSB solution is self-consistent if the 
corresponding energy exponent is positive, 9 — Op — 
d—2-\- 2C,F > 0. This occurs in dimensions 4 > d > 2 for 
both models and any 7, as well as in dimensions d < 2 for 
sufficiently long-ranged disorder in model II [7 < 7c (d) = 
2/(2 — d)], including model / in d = 2. 

The replica symmetry is broken for small IR cut-off 
TO < TO-c, where tUc = mc{T) is the temperature depen- 
dent Larkin mass, which is determined by the instability 
of the replica-symmetric (RS) solution: 

1 = 4B"(2r/i(TOe))/2(TOc), (17) 

and decreases as a function of T from rndT = 0) = 0(1) 
to zero as r ^ 00. Hence, in that case, an unconstrained 
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FIG. 1: Phase diagram for model I in d = 3 for A/rric S> 1. 
The phase transition as given by (|17|l is everywhere continu- 
ous, and the glass phase exhibits continuous RSB. The tran- 
sition temperature Tc diverges as m ^ 0. An unconstrained 
system (m = 0) is therefore always glassy. A similar phase 
diagram applies to d > 2 for both models I and II, and 
to d < 2 for model II with sufficiently long-ranged correla- 
tors, 7 < 7c (the temperature scale being set, for d < 2, by 
/i(mc(r = 0))). 



system (m = 0) is always glassy, while a strong confine- 
ment TO > wiciT) leads to an ergodic (replica symmet- 
ric) high temperature phase. This is illustrated in Fig. [T] 
where we plot the phase diagram for model I in d = 3. 
Since the temperature always enters in the combination 
2r/i, where /i is dominated by the UV cut-off for d > 2, 
we have introduced the rescaled temperature: 



f = 2T/i(to = 0) 



4TA 



d-2 



(47r)'^/2(d-2)r(d/2)' 



(19) 



for a circular UV cutoff in k space. 

The selfenergy function cr(u) of a continuous RSB so- 
lution generally interpolates continuously between two 
plateaux at small and large u. For both models (I) and 
(II) (with 7 < 7c in d < 2) cr(u) takes the form§^: 

\2-e 




a{u) = 



Here rric = rac{T), and^ 

Uc = ATml 
u„i = ATm 

where, with /„ defined in (fTSl) 



u < u,„, 

U„j < u < Uc 
Uc < u < 1, 



A 



hirric 



AA, 



■h{mc) 



(model I), 



(20) 
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(23) 



4Ad 



^ = ^(l + ^)(i^) > (modem), (24) 
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FIG. 2: Phase diagram for model I in d = 2 (with a UV cut- 
off A = 10). The phase transition is everywhere continuous, 
and the glass phase exhibits continuous one-step RSB. The 
transition temperature has a finite limit Tc —> 2Tr as m —> 0. 
A similar phase diagram applies to model II with critical long- 
range disorder 7 = 7c in d < 2. The transition line has a 



cusp-like behavior as m — > 0, 
(and ~ l/log(l/m) in d = 2) 



[Tc(0) - Tc(m)]/rc(0) ~ m^ 



and Ad := e/j^(l + P)-^ = 2r(3 - d/2)/(47r)^/2. 

In the sequel the value of cr(0) will play a central role, 
and we give an explicit expression for later reference: 



cr(0) = cr(u„i) = 



9 AT 



(25) 



In the case of sufficiently short-ranged disorder in d < 2 
(model I or model II with 7 > jcid)) the glass phase is 
described by a 1-step RSB solution which is fully char- 
acterized by three numbers: the break point Uc, and the 
self-energy parameters 



ctq = ct(u < Uc) =cr(0), 
0-1 EE cr(u > Uc) = cr(l). 



(26) 



The borderline between 1-step and full RSB is charac- 
terized hy 9 = 6f = 0- 7 = 7c for model II in d < 2, 
or model I in d = 2. In this case the 1-step solution 
can equivalcntly be obtained as the limit of a continu- 
ous RSB solution, which entails that the 1-step scheme 
is only marginally stable. In this limiting case of con- 
tinuous RSB, TOc(T) remains a decreasing function of 
T, However, it vanishes at a finite T = Tc signalling a 
continuous glass transition for the unconstrained system 
(m — 0) at Tc- The phase diagram for the marginal case 
of SR disorder in d = 2 is shown in Fig. [2l 

Away from the borderline in the (d, 7) plane, the 1-step 
solution is genuinely stable. This includes in particular 
the case of the directed polymer and the KPZ problem 
(d = 1). The phase diagram in the (T,m) plane is more 
complicated: (cf.. Fig. 2] on page [^0)1 and will be dis- 
cussed together with the 1-step solution within the GVM 
in Section lIlID3l The essential difference with the cases 



discussed above is the nature of the temperature-driven 
glass transition at smaU mass which becomes a discontin- 
uous random first-order transition. Similarly to ordinary 
first-order transitions, two locally stable solutions coex- 
ist at least close enough to the transition: Here these 
are the RS solution and the 1-step solution. The rough- 
ness exponent of all 1-step solutions equals the thermal 
exponent, C = Cth = (2 — d)/2. 



2. Effective action and functional RG 

In an effort to connect the results described above to 
the functional RG approach, two of us performed a cal- 
culation of the effective action T[u] for the model ([4]), cf. 
Refs. IlllflSl One starts by defining the standard gener- 
ating functional W[j] for connected correlations in the 
replica theory with, in general, j = {j^{x), . . . , j"(a;)): 



[a](u) + wi^ f 



,vv[j] 



Z[j]^l[Zv[n 




FIG. 3: The MP function [a]{u) 



for m > 0. Changing 



m , only the lower plateau will move, while the remainder of 
the function (for u > u,n) remains unchanged. 



(27) 



Zvif] = / PMe-'3^^'["l+■/■=^^'°•"^ 



(28) 



and the effective action functional is defined via the Leg- 
endre transform: 






5W 



[j] 



(29) 



In Refs. ITlllTsI the effective action for a uniform configu- 
ration r(u) = ■^r[{u'^ = u""}] was studied in the large- iV 
limit. It has a non-trivial limit in the regime 

u^ r^ N, r ~ iV. (30) 

It was further computed in an expansion in the "number 
of replica sums" (which is effectively a cumulant expan- 
sion) as follows: 



nu) = E 



2T 



^E^("") 



2r2 



^.T^s^'Hu^"! 



6r3 



(31) 



ibc 



ab 



up to a constant. Here and below u' 
^abc _ yO^yfc^yC^ This defines unambiguously (for any 
N) the second cumulant of the renormalized disorder po- 
tential R{u) = Rm{u), where we will usually keep the 
TO-dependence implicit. At large N, it obeys the scaling 
form 



R{u) = NB{u^/N). 



(32) 



Remarkably, the scaling function B{z) satisfies a closed 
equation: 



B\z) = B' i 



B' (z + 2Th + 4/2 



B'{z)-B'{0) 



(33) 



where z — vF'/N, and /„ — Inirn) was given in p8|) . 
As we will show later, this equation is only valid in the 
non-glassy regime since it was derived under the im- 
plicit assumption that the replica symmetry is not spon- 
taneously broken as u — *■ 0. Taking a derivative with 
respect to the mass m and introducing a scaled func- 
tion b{x) = 4Adm'^''~'^ B{xm~'^'') a FRG equation was 
derived for the renormalized and scaled potential corre- 
lator b{x) (denoted b{x) in Refs. [lllllSl ). Remarkably, 
this FRG equation admits a natural continuation to the 
glassy regime. In that continuation a linear cusp exists 
for all m < rric. The resulting large- A'^ equation, valid 
to all orders in e = 4 — d, exactly matches the equation 
previously derived to first order in e (1 loop) but for ar- 
bitrary A^ by Balents and Fisher^ by a quite different 
method. In addition, for all universality classes charac- 
terized by ((i,7), such that 9 > (ensuring continuous 
RSB in the GVM), the fixed points b*{x) could be re- 
lated to the solution of the GVM. All the above strongly 
suggested that the derived fiow equation was the correct 
large- A^ limit of the FRG. 



3. Relation between GVM and FRG in the case of 
continuous RSB 



However, establishing a link between the GVM results 
and the FRG is rather subtle. To a considerable extent, 
this has been achieved for the case of continuous RSB in 
Refs. [IIII13I , while the same task for situations exhibit- 
ing 1-step RSB has remained an open problem, which is 
addressed in the present paper. In order to review the re- 
sults for the continuous case, we recall the multiplication 



formulae for two ultrametric matrices AB = C: 

A^B" = C^ (34) 

{A^ - [A][n)){B^ [B][u)) = C^ - [q(u), (35) 

A{Q)B'' + S(0)A" = C(0), (36) 

where: 



A" -^ A 



aa I A{u)du, 





[A]{u) := uA{u)~ / A{u')du' 
Jo 



(37) 
(38) 



We also introduce the notation for the diagonal clement 

A := Aaa- 

From (|20|) above one finds: 

{m/mcf, u < u™, 
m^ + [cr](u) = m^ X <( (u/uc)^, u„ < u < u^, (39) 

1, u„ < u < 1, 



and from (|14p we obtain 



G=(fc)-[G](fc,u) = — ^ 



T 



(40) 
(41) 



.9.T^ + M(u)' 
G{k, 0) = G(fc, u„) = Ta{u„,)gl (42) 

Further, using ^^(u) — u^cr(u), we find: 

(7(u) du 



G{k,u)=Ta{um)9i+T 



u,„kT^ + M(u)F 



(43) 



(we often denote & = ^a here and below) which relates 



-1 - 1.2 



G{k, u) to the self-energy ct(u), which for g^, = k 



is given by ^ and ^ 

In the GVM the 2-point correlation function is thus 
given by contributions from states at all distances u as: 



G(fc) 



gk + cr(u™)5^ +gk 



^ du [cr]{u) 



(44) 



Using (|39p and integrating over u yields the final result: 
G{k) - Cm-^'^+^^^f{k/m), (45) 

in the limit to -> 0, with /(O) = 1 and f{z) ^ ^^(''+20 
with C = C-Fj exhibiting the anticipated scaling. 

On the other hand the FRG allows one to compute the 
two point function using ([3T|l through the exact relation: 

(27:)' S^'Hk + k')G-J{k) = f^\u=o , (46) 

which, upon matrix inversion yields the zero-momentum 
correlation function: 



G,. = o, = -H|M,i,. 



m* m^ 



(47) 



It turns out that this does not yield the full MP result 
(|44|45|) . Indeed, the contribution from the integral over u 
in (|44|) is missing. However, Eq. (|47|l exactly reproduces 
the first two terms in (|44p . i.e., the contribution from 
the most distant states corresponding to < u < u,„, as 
well as the trivial term T /m? describing the connected 
correlations G'^{k — 0). Using the natural continuation 
to the glassy regime of the RG fiow associated to the 
solution of (1331). it was indeed inferred that 



-2B'(0) = To-(u^ 



(48) 



as will be derived rigorously in Section IIVI This is al- 
ready a non-trivial result: obtaining it within the GVM 
does require the RSB Ansatz, while the FRG considered 
in Refs. [lllllSl does not make reference to ultrametric 
matrices and their properties. 

The fact that the FRG calculation of Refs. [Hil re- 
produced only the contribution of the most distant states 
within the Parisi hierarchy was argued to arise because 
r and W were computed by introducing sources which 
split all replicas, i.e., u"'' ^ (and of order y/N) for all 
a j^ b. However, it remained to be understood how the 
full result (|15|) can be recovered and how the crossover 
to spontaneous RSB takes place as u"'' -^ 0. 

This is the aim of the present paper. We will show 
the existence of a regime whose large- A^ limit is different 
from dSOl): 



N°L 



-d 



r~ ATU 



(49) 



and where RSB cannot be neglected. This will allow us to 
recover the complete GVM result, i.e., to correctly per- 
form the (non-trivial) zero-field limit u^ -^ in ([46]) . We 
first show, in the following section, that it is equivalent, 
but physically more transparent, to study the generating 
functional W rather than F. Further, we will revisit the 
attempt by Balents, Bouchaud, and Mezard^ to connect 
RSB and FRG. Although we choose a different observ- 
able than those authors (so as to connect more directly 
to the FRG) we employ techniques similar to theirs, at 
least in the regime v^ = 0(1). 



C. Observable studied here 

The observable studied here is the free energy V\v\ of 
the elastic manifold in a quadratic potential well centered 
at position v^, and defined by2£: 



e-^M/^= /l?Me-^-["^^l/^, (50) 

Hv[u;v] = 2 / 3fc Vfc -'^fcP + / V{x,Ux) , 

with gk = k^ + TO,2. For a uniform v^ = v, Vk ^ 
v{2T:Y5^'^\k) onchasSl: 

Hv[u]v\ = Hv[u] ~ m^v u^ + -L'^rn^v^. (51) 

Jx ^ 
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It is then clear from the definitions (|27l50p that the statis- 
tics of V[v] can be obtained from the functional W[j] 
with: 



Jk 



9k '"k 
T 



(52) 



Hence we denote W[v] — W[j = {g ^v)/T], such that: 



JV[v] = g 2T A Sfc ' S:,. i-Mfc TT (,-Vlv^]/T 



(53) 



It is thus clear that the functional can be expanded in a 
cumulant expansion, i.e., in the number of replica sums: 

w[v] - ^ E / Sk'vivi, + ^ E ^[^"'i 



s^E^^^^[-"i 



6T3 



(54) 



a be 



up to a constant, where again w"^ — v'^ — v^, etc. Each 
term in the expansion is a functional, i?[ti''''] = i?[{ti°''}] 
with fixed a, b, etc., and represents the free-energy cumu- 
lant s: 



V[v]V[v'] =R[v-v'], 



(55) 



and similarly for higher cumulants (the overbar denoting 
the disorder average over the random potential V). For 
a uniform configuration Vx = v one has 



R[v] = L'^R{v) 



(56) 



which defines R{v). Note that with arguments in brack- 
ets [...] we denote functionals, while (...) is reserved for 
functions. When discussing uniform v we usually switch 
to the function i?(w), separating out the volume factor. 

As shown in Refs. [48|[49l by performing the Legendre 
transform (|29p . this observable is directly related to the 
function R{u) of the FRG: 



Riv) = Riv), 



(57) 



i.e., the two functions are the same, and this holds for the 
functionals, too. Therefore, by computing R[v], which 
we do here for large TV, we simultaneously compute R[u\ 
as defined from the effective action (|3T|) . The fact that 
([55)1 defines an observable which is easy to measure has 
allowed for a numerical determination^'^ of R{u) or, more 
precisely, of the force correlator A{u) — —R"{u), for N = 
1 interfaces at T = 0. 

The derivative of R{u) at T — contains information 
about the shocks. Indeed, if one computes the ground 
state u{x; v) for a fixed well position and defines the cen- 
ter of mass displacement u{v) := L~'^ J u{x; v), one finds 
that the latter exhibits jumps as v is varied. The statis- 
tics of these jumps is encoded in the functions R, S etc., 
for instance one has: 



[0{v) - v^[ui{v') - v'^ = m-'^L-'^Aij{v - v'), (58) 



and the cusp of Ay (v) = —didjR{v), i.e. the derivative of 
A at argument 0, is proportional to the second cumulant 
of the jump sizes. 

We now turn to the calculation of this observable in 
the large- A^ limit. We will perform separate calculations 
for the two scaling regimes v'^ ^ 0(1) and w^ ~ 0{N) in 
the next Sections and check that they match. 



III. REGIME OF SMALL v, vl = 0(1) 

We want to compute the generating functional ((27|) in 
the form 



W\v\ 



n^ 



v\]'' 



3k — 



9l^< 
T 



(59) 



This can be carried out at large N through the saddle- 
point method as in Ref. l46l introducing the term 
iXabiNxab — Ua ' i*b) and integrating over the field u. This 
is detailed in the next Section where we study the case 
where w^^ ~ 0{N) which distorts the saddle point away 
from (|14ll5p . In this section, wc study w^^ = 0(1) and 
the above saddle points are unchanged (more precisely, 
they are shifted only by terms at most of order 0(1/A) 
which are discarded). Since there is spontaneous RSB, 
there are in fact many saddle points equivalent under 
replica permutations. Hence the above average must be 
written as a sum over all equivalent saddle points: 



n^v'b-'^] = C„E(e 



T.a L 3° 



(60) 



where n € Sn belongs to the group of permutation of n 
replica which is used to label the saddle points. Around 
each saddle point the measure is Gaussian with correlator 



«u: 



Glbix-y) := G^(Q)^(6)(a;-y) where Gabix) is 



the ultrametric matrix given by (I14p . Hence we obtain: 



^W[v]~W[0] 



E. 



exp 



1 

2r2 



T. 9k'G:b{k)vl,-vl 



(61) 

where ^^ denotes a normalized average over permuta- 
tions, i.e., X^Tfl ~ 1- T^^^ remainder of this Section is 
devoted to the analysis of this formula. 

In principle, by expanding this formula in powers of v 
to all orders, and regrouping terms one could check that 
it is indeed possible to put it into the form ([51]) . and to 
compute all the cumulants (we denote by S*^"-* the n-th 
cumulant). This is a formidable task, however, and we 
will focus here only on the second cumulant function R. 
Before computing it directly, let us give a fiavor of the 
direct expansion in powers of v. 



A. Direct expansion in powers of v 

The expansion of (|6ip in powers of v starts as: 



W\v 



Z ^-j^ 



—y 

ab 



9l^GUk)v% ■ vl 



(62) 



Let us recall the structure and parametrization of a hi- 
erarchical Parisi matrix. Dropping temporarily the k- 
dependence one haa^ii^^ 

Gab = [G- G{n)]5ab + G(n)(]l„)a6 , [RS], (63) 
Gab = [G - G{l)]5ab + [G(l) - G(n)](l(';))afc 

+G(n)(l„)afc, [1-stepRSB], (64) 

dG(u) 



Gab = [G- G{l)]5ab 

+G(n)(l„)„fc, 



d\i ■ 



du 



<^[r\^ 



[continuous RSB]. 



(65) 



We have defined ]l„j as the matrix made of n/m identical 
blocks along the diagonal, each block being the mhy m 
matrix with all entries equal to one, and !„ ^ = !„. As 
an example, 

/111000000000\ 
111000000000 
111000000000 
000111000000 
000111000000 
„(i2) _ 000111000000 
000000111000 
000000111000 
000000111000 
000000000111 
000000000111 

yoooooooooiii/ 

With these formulae one easily checks that 

G"" := y^Gab = G- I G{u)du + nG(7i) 



(66) 



n^O ^ 



G(u)rfu. 



(67) 



For symmetry reasons, the permutation average must 
yield a replica symmetric matrix 



J2 G:,^aSab + p. 

' ^ TV 



(68) 



Setting a^h one finds /? = (n — 1) ^ X^b^ti Gib — {G'^ — 
G)/(n— 1), while multiplying with 5ab shows that a + /3 = 
G. Thus, one finds: 



^.^-^ " T 



1 



QC 



(G" - nG)5ab + G-G'' 
n(G - G^)l 5afa + G - G^ (69) 



which is a replica symmetric matrix by construction, but 
with non-trivial entries in the case of RSB. Terms of 
higher order in n have been neglected, except for the 
linear term in the replica diagonal part which we retain 
for later use. 

Inserting in (|62p we finally have 



W[v] - W[Q\ 
_ 1 
~ 2T2 



(70) 



E 

a 

—y 



9-.' 



a\2 



2T 



ab 



G"(fc)+n(G(fc)-G"(fc)) 
g^'iGik)-G^ik)K-v'_, + 0{v^) 



On the other hand, expanding (|54|) we obtain, for n = 0: 

1 



W[v] - W[0] 



-y[(- 

+0{v^), 



b 

Tg~k 



T2 



Rm]vk-vt, 



(71) 



since the fourth cumulant must start as w^ and the third 
as v^ (as one cannot construct STS invariant combina- 
tions of smaller degree - see Appendix B of Ref. 22. for an 
argument in the case A^ == 1). Here and below we adopt 
the notation 



Rrr 



g»fcx+zfc'x'^»^^j^j 



S^Rjv] 
i2TrfS{k + k')R'l[v] 



(72) 
(73) 



where the last line uses translational invariance and one 
has set v^ — v after taking the derivatives. Further, 
we have assumed 0{N) symmetry so that the second 
derivative of R is diagonal and independent of the 0{N) 
index i in (|72p (where no index summation is assumed). 
Identifying ^Q and ^1^ shows that for n = 0: 

Gik)^Tgk-glRm , G=(fc) = Tg^, (74) 

the second identity being a simple consequence of the 
STS symmetry. Hence, once the functional R = R oi the 
FRG is known, the correlation function computed via the 
GVM can be retrieved as: 



J2 Gl,{k)^Tgu6ab-glRm- (75) 



-k^-k/TT 



This property of R, here guaranteed by its definition (pT|) 
via the effective action, does not hold for the observable 
defined in Ref. |47| (having a similar form but with G -^ 
G~^). This made the comparison of their results with 
the FRG problematic. 

To evaluate the identity ((7i|) at A; = 0, we use that for 
a uniform v, d/dv = J dxj:^, and thus: 



R'LM = ^, 



R!.L\v\=R"{v), 



(76) 
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which yields: 



G(fc 



T 

= — 



R"{0) 



(77) 



This relation is exact, and the task is hence to evaluate 
R"{0). There is however a crucial subtlety in evaluating 
this derivative at t; ^ 0. If one uses that R{v) — R{v), 
together with the result of Refs. [ll|[l3l : 



R{v) = NB{v^/N) 



N 



(78) 



a relation also obtained directly in Section llV Bi one finds 
that ([77|) coincides with formula (|47|) . However, this re- 
sult is valid only in the region v^ ^ N and, as pointed 
out above, it does not reproduce the full MP result for 
G. To obtain the latter, as we show below, one needs to 
be more careful in the i; ^ limit and compute R"{0) 
in the region where v"^ ~ L~'^ ^ N. One could say that 
the u — !■ and N —^ oo limit do not commute, or more 
accurately, that to obtain contributions of all ultrametric 
states (and not just the most distant ones) one must take 
the limit u — > with great care. 

The next order 0{v'^) is obtained in App. [B] by the 
direct expansion method and in the next Section by a 
more powerful method which can handle all orders, and 
to which we now turn. 



B. Second cumulant from a tw^o-group analysis 

While the results of the previous section were com- 
pletely general and independent of the RSB scheme, we 
now focus on a specific choice of external sources. In or- 
der to compute the second cumulant function, one best 
uses two sets of replica, which we denote by v^ = v^, 
for a = 1, . . . , n/2 and w" = u^ for a = 1 -|- n/2, . . . , n. 
Inserting into (|54p we find 



WW 



where 



-Wm^^^Jj^'{vl-v'^, + vl-v'_,) (79) 



1 n' 
'2T2^ 



{R[v^^]-R[0]) + O{n^), 



,,21 



(80) 



since all higher cumulants yield higher powers of n. We 
compare this with expression (|6ip . slightly rewritten as: 



^wiv]-w[a] 



V exp 

^ ^TT 



(81) 



ab 



W^J2 9k^Gab{k)v 



k ' ^fc 



In order to perform the sum ov er p ermutations, we in- 
troduce Ising spins t" as in Ref. |43: 



,,21 



(82) 



where t° = —1 if 7r(a) G [l,n/2], r" = +1 otherwise. 
Each configuration {ra} is left invariant by [(?i/2)!]^ per- 
mutations, the number of distinct spin configurations be- 
ing Cn = "■■/[("■/2)!]^- They all correspond to vanish- 
ing total magnetization, i.e., J2a=i '''"' — 0- Inserting ([82]) 
into ([5T|) the term linear in r vanishes, and one finds: 



,WM-Wm ^ exp ^-^y^.g-2GC(fc)(|,l|2 + |,2p) 
*-"" {t} L '^ ab 



Gab '■— Gab ~ G Oab 



(83) 



The first factor equals the term proportional to n in (j79p , 
as seen upon using (j74p . The prime on the sum indicates 
that the sum extends over all Ising spin configurations 
subject to the global constraint J22=i '''^ ~ ^- Identifying 
the above expression with (|79p . we arrive at the formula: 



R[v] - R[Q] = hm 



4T^ 1 






,n/2 



(84) 



xE 

{r} L 



exp 



/>i-i' 



^G„j(t)TV' -1 



Note that we have simplified notations by renaming 
v'^^ — > V. At this stage we can check the small-u ex- 
pansion again: 



RmH\ 



(85) 



*o2n2 



9k 



li-:^ i^klO .. G.,(fc) 



2Y^ 



^a—h 



2T2 



^^{r} 



If we denote the normalized average Aab '■= "^i^xTaTb — 
„^2 ^r 1 TaTb then one has Aau — 1- Further, the iden- 
tity ^ J2b ^Q& = ^aa + (n - l)Aa^b implies Aa^b = 
1/(1 — n) and one can thus write Aab = (^ — nSab) / [^ — n) 
and EafcGab^afc = -ny{l-n)[Gik)-G^ik)]. Using this 
in ((85|) we recover the result ([74]) of the previous Section. 
Note that the same result is obtained from the piece ^ n^ 
in the replica-diagonal part of (|7D)) , while the two-replica 
sum vanishes. 



To evaluate the restricted spin sum in ([M| we use the 
same method as in Ref. [53, leading to Parisi's nonlinear 
diffusion equation in the form discussed by Duplantier'"'^.. 
We first eliminate the constraint of zero magnetization 
and rewrite (IMl) as: 
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R[v] - R[0] = lim 



4^2 2»+ir(-n) 



n^o n2 [r(-n/2)]2 



dy 



{r} 



exp 



— 2i 12 '^ 

5fc kfcr Y- 



8T2 



a6=l 







(86) 



which makes use of an identity derived in Ref. |59J, vahd 
for n < only. 

The evaluation of a spin sum such as it appears in (I86p 
is standard in the mean-field theory of spin glasses. Here 
we summarize the main steps following Duplantier— . We 
write the matrix which couples the spins as 



1 



with g"^ : 
q(u) = 



Qab = -^ / 5fe \Vk\ Gab{k), 

J2b^ab ~ 0. Using (|i5|) this becomes 

a{u)du 



(87) 



4T 



\Vk\ 



o-(u„0 + 



u„ [5-i + H(u)]' 



Let us now assume that qab has a iiT-step ultrametric 
structure with breakpoints at 

n=uo ^ ui ^ ■■ ■ ^ UK ^ UK+i = 1, (89) 

and entries parametrized by q{u) with 

q{u) = qi, for u<> =:^ u -( u^+i. (90) 

Further we define q_i = 0. Let us introduce the "partial 
partition sums" : 



gi{y) = e^^'^^(^) 



H exp 

{r"},a— 1,..,U£ 



(91) 



I Er'^(ga6-?£-i)r'' + f]A 



a, 6=1 



which defines geiy) and ipdy) for £ = 0, ... X + 1. Obvi- 
ously, gK+i{y) = 2cosh(y) exp[5 - g(l)]. 

The expression (1861) can then be rewritten as: 



R[v] - R[Q] = hm 



-2T^ 



rnO n 
2^2 



dy[5o(y)-(2coshy)"] 

) 

dy[V'o(y)-ln(2coshy)] , (92) 



r(-n/2)] Vr(-n) = -4/n + O(n0). 

58l . the functions gi{y) satisfy a re- 
cursion relation, the idea being as follows. Consider one 
of the n/u£ equivalent groups of size u^, for instance the 



where we used that 
As shown in Ref. 



l,...,Uf. It contains -^^^- subgroups of 



first one a — ^ "• t+ ^.r^nfcino _Hi 

— The only coupling between the subgroups 

92 



size u^+i- 

is through the matrix Ag^l" with Aqg — qe ~ qe-v 

which couples uniformly all spins in the group. They can 



be decoupled by a Hubbard-Stratonovich transformation, 
adding a term z^/Aqg J^ILi '''a, where 2; is a Gaussian 
random field of variance Aq^, acting on all spins in the 
group. This allows to perform configurational sums inde- 
pendently within each subgroup and yields the recursion 
relation: 



9iiy) 



(9i+i\y + zy/Kqiy'^^\ , 



(93) 



where here and below {. . .) ^ '.— J_ 
the average over a unit Gaussian. 



00 dz 
2-K 



/^ denotes 



C. Continuous RSB: Second cumulant from an 
evolution equation 

In the case of continuous RSB, we have to take the 
continuum limit if —> 00 in the above: Au^ = u^+i — u^ — > 
0, q{\ii) = qi ^ q{u) and Aqg — > 0, while Aqi/Aui -^ 
dq{u)/du. For the function ?/;(u^,j/) = V^(y) -^ ip{u,y), 
Eq. ((93|) yields a differential equation: 



a,^ = -iM!)[a2^ + u(a,V')^ 



2 du 
with initial condition: 

V'(uc,y) = ^(l,y) = 



(94) 



9-9(1) 



ln[2cosh(y)], (95) 



where we assume q{u) to be constant for Uc < u < 1, as 
is generally the case, cf., (^01) . 

Eq. (IM)) must be integrated from Uc to Um, where q{u) 
reaches its lower plateau go = g(um). If go 7^ 0, the 
recursion ([93| for £ = shows that we have to perform a 
last convolution to obtain 



V'(0,y) = (i^(u,„,2/ + zv/g(0))>. 



(96) 



The solution of 
at large \y\ like: 



\y\ — »oo 

V'(u,2;) « \y\ 



with initial condition ([55]) behaves 



1 



du'g(u') 



-ug(u) 



(97) 
as can be seen by substituting ([57]) into ([M|) . This is 
also confirmed by noting that for \y\ 3> 1 the sum in 
([86|) is dominated by the configuration with all t'^ — 
sign(y), which leads to the simple approximation ((97|) 
for V'(u,y). This behavior implies •il}{um,y) ~ \y\ + q'^ — 
|j/|, since g^ = from the condition G"^ = 0. Hence 
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the integral ([92)1 converges, the subleading terms of ip 
decaying exponentially at large \y\. 

We can now state the main result of this section. In 
the regime v^ = 0(1), the second cumulant functional 
can be obtained from the saddle point of the GVM as 
follows: 

/oo 
dyy[M{Q,y)-t'^nh{y)] ,(98) 
-oo 

where the function M{u, y) = dy'ip{u, y) is the solution of 
duM = _iM!^ (a2M + 2uMdyM) , (99) 



2 d\i 

M{uc,y) = tanh(y), 

in the interval u G [u„i, Uc], and 

A/(0, y) = (m(u™, y + z^/^) 



(100) 



(101) 



Eq. (|98p follows from Eq. ([5^ by integration by parts, us- 
ing that limy^ioo^ {V'lO, y) — ln[2cosh(y)]} = as dis- 
cussed above. 

From the above and (IST)) we see that the dependence of 
the functional R[v\ on the field Vx occurs only through the 
combination q{\i) = -^^ J h{x — y, u)v{x) ■ v{y), where 

h{k,u) — g^ G(fc,u). Hence for a uniform Vx — v one 
easily sees that q(u) ^ U^v'^ and we can expect that ([98]) 
will assume a scaling form with scaling variable vL'^^'^. 
This dependence on the system size is very different from 
the one in the regime v^ = 0{N) and will be commented 
on further below. 

The convolution (jlOip only results in an additive con- 
tribution to the potential correlator. This can be seen 
using the identity: 



dyy 



-^[$(oo)-$(-oo)]. 



(102) 



proven in App. \C\ for functions <i> whose derivative de- 
creases faster than l/|y| at infinity. Applying p02p to 
<I>(j/) = M{um,y), we can rewrite (jM)) as: 



R[v]-R[0] = -^J g^^G{k,u^)\vk\' 



-2T 



(103) 

pOQ 

2 dyy[M{u,n,y)-ta.n-Hy)]. 



This result can also be derived via an alternative route, 
providing a useful check of the above formalism: In 
(|86p we can directly separate out the contribution form 
the most distant states by writing Gab{k) = Gab{k) + 
G(fc,0)(l - ndab), with G\k) = 0, and G(fc,0) = 0. 
One can easily check that the piece G(A:,0)(1 — nSab) 
produces the first term in (I103|) . The remaining part 
leads to the same expression as ([5S)) . with the replace- 
ment G —> G. The only difference in the subsequent 



evaluation is that there is no need for a convolution in 
the end, since G^ik) = 0. This estabhshes (fT03l) . 

Note that the first term in (jlOSp has the expected form 
for the contribution from the plateau < u < u,„, cf., 
Eqs. (|44l48p . with a coefficient of \vk\'^ 



- -;^9k'^G{k,\i„,) 



1 



ra(u„0 - B'(0), (104) 



independent of k. Clearly, this term is the only contri- 
bution in the case of a replica symmetric solution. We 
now show that it also gives the dominant contribution at 
large v"^ of order 0(iV"), but u^ > L^'^. 

1. Limit of L'^v^ :^ I 

In this case, Eq. ([M]) can be rewritten as 

1 . // / N dM 

- - (M" + 2uM'M) = » (105) 

2 dq 

since q is large (from now on we denote dy by a prime). 
One can integrate this equation to M' + u{M^ — 1) = 
0, where the integration constant (w.r.t. y) is fixed by 
the fact that M{uc,y) -^ 1 for large y, which cannot be 
changed by the evolution ((99)) . The solution of (|105p is 
then Af (u, y) — tanh(uy). One can check that the flow of 
Af is attracted to this simple "fixed point" as u — > u„i, if 
w 3> w*, with v^, defined below. In that case, the second 
term in (|103p becomes: 

/oo 
dy y [tanh{umy) - tanh(y)] 
-oo 

-11/ dyy[tanh{y)-l], (106) 



4T^ 



and the large-w behavior of the second cumulant in this 



0(1) regime is thus: 



R[0] - R[v 



ti^oo Tcr(Um) 



71" ^9 / 1 



(107) 
The leading behavior is quadratic in v and corresponds to 
the contribution of the most distant states to the full (in- 
verse) correlation. Hence it should match the result ob- 
tained in the FRG for the regime v^ ~ N, cf. Eqs. (|47l48p 
or equivalently, Eqs. (|77l78p . Indeed, it does! 

For a uniform Vx = v and in the limit T —^ 0, one finds 
for both models dHEl), using (|20l22p : 



R[0] - R[v] 



^2-e^2^d ^2 



-26» 



(108) 



This suggests a crossover around v '^ v^, defined as: 

m^L'^vl = m"^ ^ v^ = [mLY'^l'^mr^ . (109) 

The physical meaning of this crossover scale will be 
discussed in more detail in the context of the models 
with one-step RSB below, cf. Sec. IIIIEI 
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2. Perturbation expansion for L v^ -^ 1 

Since for g(u) — 0, one has M{u,y) = tanh(y) =: 
mo{y) there is a uniquely defined expansion in powers 
of q, which since q r^ v^ is equivalent to the direct expan- 
sion in powers of v'^ of Sec. IIII Al 



M{u,y) = mo{y) + mi(u, y) + m2(u, y) 



(110) 



Each mn{u,y) contains only terms of degree g". They 
satisfy the recursion: 



1 



m„ 



^(u) 



/n-l 



Ci + u E^ 



^p''^n—p—l 



\p^0 



(111) 



where here and below dots denote du- The initial condi- 
tions are TO„(uc,y) = 0. The final result for R[v] is: 

/C30 
dyymr,{u„„y) 

= :J2Rn[v], (112) 



where we recall that 4T'^q{u) = /^.^^ ^|wfcpG(fc, u). 
Hence Rn contains all terms of degree q" ~ w^", and 
we are thus effectively computing the derivatives at the 
origin, 7?(2")[„ = o]. 

This calculation is performed in Appendix [D] and the 
results are indeed consistent with the direct expansion of 
Sec. IIII Al although the method is quite different. The 
lowest order term reads: 



RiM 



-2T^ / q{u)du 
Jo 
1 



,-2 



o 9k' 



G{k) - G"(fc) 



wfe 



(113) 



using that /^ G(u) — G—C^, recovering the result 
which yields, at small v, the full result of Mezard-Parisi 
for the correlation function. 
The next-order term is: 



R2\v] = -T' 



du q (u) — 




(114) 



As discussed in Refs. |8||22||38||48!|49| the fourth deriva- 
tive at zero of the FRG function i?'^'*^ [0] is a direct mea- 
sure of susceptibility fluctuations. Indeed we find that 
our present result has the general form of susceptibility 
fluctuations within a Parisi Ansatz, defined and derived 
in Ref. \M ■ 



3. Thermal boundary layer: General formula 

It is possible to resum the derivatives, order by order 
in temperature and derive the thermal boundary layer 



(TBL) form: 
R[v\ 



Rn - -2 



5fe^G'T=o(fc)|wfe| 



+T^fi[v]+T^f2[v] 



(115) 



where the first term is the "zero-temperature limit" of the 
leading small-ii quadratic term^^, and the higher-order 
terms in the expansion in T are scaling functions of the 
boundary-layer variable: 



Vk 



Vk/T, 
q{u)^\j^gl^\vk\^G{k,u). 



(116) 
(117) 



This structure, which appears already in the 1-loop FRG, 
can be computed exactly in the large- A'^ limit here. 

This structure appears already in the 1-loop FRG 
where one finds that at T > the non-analyticity of 
R[v\ is thermally rounded i^^'^^i'^^ Since temperature is 
irrelevant in the RG sense when 9 > Q this rounding oc- 
curs only in a layer determined by wm^ ~ Tm^ , which 
becomes smaller and smaller as T ^ 0, hence the name 
thermal boundary layer. Here, for the first time we com- 
pute its exact expression in the large- A^ limit. 

We now perform a (slightly more formal) expansion 
which allows to obtain (|115p . The idea, looking at (|99p . 
is that the nonlinear term contains an extra factor of u 
and that for u^ < u < Uc one has u ~ T. Hence it is 
natural to expand in the nonlinearity to generate a low- 
temperature expansion. Of course this must be checked a 
posteriori. Further, we find it more convenient to use q{u) 
instead of u to parameterize the ultrametric distance. 
Thus rewriting ^ as dqM = -\ \M" + \x(M'^)'\ , and 
expanding M = Mq + Mi + . . . , formally in powers of u , 
we have to solve the hierarchy: 



dqMo 

dqMl 



a,M„ = 



- M" 

1 



- (Aff + 2uMoM^) , 

1 / ""^ 

-- M;' + 2uEMpA<_i_ 

\ p=0 



(118) 
(119) 

(120) 



Note that this hierarchy is formally similar to (|llip . 
apart from a shift in the index of the linear term on the 
RHS which indicates that we actually perform a different 
resummation. 

The initial condition is AIo(uc,y) — tanh(j/). 



Af„>i(uc,y) = 0. 
fine: 



Again, primes stand for dy. We de- 



R[v] -i?[0] = na[v]+ni[v 
7^oH = -2T\{0) 



(121) 
(122) 



7^„>lH = 2T' 



dyy[Mo{q{Um),y) - tanh(y)] , 
dyyMn{q{um),y), (123) 
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where we use the notation TZn to distinguish from the ex- 
pansion in powers of q discussed in the previous Section, 
and different from the r^"^] defined in (jllSp . 

The equation for Mq is a simple diffusion equation with 
solution: 



Mo{q{u),y) = (tanh(y+V'?(uc,u)z)),, (124) 
9(Uc,u) := q{Uc)-qiu). (125) 



Thanks to the identity (|102p . when Mo{q{um),y) from 
p24p is substituted into p22p . we obtain 



7^oH = -2r2g(u,), 



(126) 



whereby the terms proportional to q(0) = q{um) cancel. 
In the limit T ^ ()126p becomes the same as the first 
term in (jllSp . since: 



G{k) = G%k) + il-Uc)Gik,Uc)+ du G(fc, u) 

^^"G{k,u,). (127) 

Here, we used G^, Uc ~ T and the fact that G(fc, u) has a 
finite limit as T — j 0. 

To obtain the next order term, TZi [v] , we solve 
Eq. (fTT9l) : 



Mi{q,y) 



qi^o) 



dq'dy'V{q,y;q',y') 
xu(g')[(MoAf^)(g',2/')] 



(128) 



where: 
'D{q,y;q',y') 



W-q) 



VMq' - q) 



exp 



[y-v') 



J\2 



2(9' - q) 



(129) 



is the (reverse) diffusion kernel satisfying {dq + ^dy)!) ~ 
S{q~ q')5{y~ y')i ^^d u{q) is the inverse function of (/(u) 
for Urn < u < Uc- Eq. (|128p can be rewritten as: 

fg(uc) 



Mi(g,2/)= / dq'u{q')y. (130) 

Jq 

(^{MoM'^){q\y + z./Y~q))^, 
which inserted into (|123p gives: 



7^l[^'] = 2r2 / dg'u(q') X 

"'9(0) 



(131) 



X / dyy{{MoM'^){q\y + z^q'-q{m. 



We can now use the identity (|102p for $ = MqMq which 
has rapidly decreasing derivatives and satisfies $(±00) — 
0. This yields: 



7^l [v] - 2r2 / dq u{q) / dy y(MoM^)(g, y) 

Ja(0) J-00 



Integrating by parts over y and using the solution (|124p 
for Mq one finally obtains: 



q(^c) 



TZi[v]=T^ I dqu{q) I dy 1 - (tanh(?/ + 2;i^g(uc) ~ q)ta,iih{y + Z2\/q{uc) - q))zi,z2 

Jq(0) 



(132) 



r 



Note that we have used (Mq — 1) as primitive of 2MoA/g, where: 
since we need it to vanish for large argument. Using 



i/'(a — 6) := / dy [1 — tanh(y + b) tanh(2/ + a)] 

= 2{a-b)coth{a-b), (133) 



m 



T 



(136) 



has a finite limit as T — > 0, where we introduced the 
we obtain the exact scaling functional for the thermal definitions: 
boundary layeB^l: 

7^lH=TVl[«], (134) 



q, := 9(uc) = 9(1) = - / gf\vk\^Gik, u,), (137) 

'J k 



ri[v\ = 2 / dqu{q) x 



I \ Ir^X ^"'(0) / I- |2 

(j™ := q{\\rn) = q[0) = — -. / \Vk\ , 



z^/2{qc ~ q) coih{z^y2{qc - q)) 



(135) 



9 - 7 / 9k'\vk\'G{k). 

•I K 



(138) 
(139) 
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4- Thermal boundary layer: Uniform v and scaling 
function for m — > 



Here and below we define: 



Let us consider tiere a uniform Vx = v for wliicli one 
has: 



9(u) = ^^^^G(fc = 0,u). 



G(u) = G(fc = 0,u = Tu) 



(140) One can rewrite ri[v] = L ri{v) with 



(141) 



m'^v^L^ 



ri\v\ 



dGu{G) 



\v\L'^/^ 



V2(Gc - G) coth 



|{>|L'^/2 



V2(Gc - G) 



(142) 



where here we denote Gm ■— G(um) and Gc := G(uc) 
(not to be confused with the connected correlator G'^). 
Let us now specify to models I and II. One finds: 



G(Q) 



1 



Ai 



A{A - 612) m2+e 2 + 
where A was given in (j23p . Hence: 

Am^(4-26l)^ 



'u^+^ 



(143) 



Q(G) 



[8-(4- 



with Uc = Am^, Am^+eG^^ 



92)ATO2+eG]2+« 

and Am^+^Gc 



(144) 



2 

2-e 



T=m-^^^) ■ Using (HH we obtain the TBL func- 
tion as a double integral in (|142p . In the limit m ^ rric 
this simplifies to: 



r-y[v] = m^v^L'^H, 
2 



1 



V 



He{x) 



1 
2 V2 



2A f *"^' 

2+7 rw dt 





^2 + 8 



(145) 
zx\Jt coth(2;a;Vr) 



where w* = L^''/2,^j-i-2 is the scale obtained in (|109p . 
Further, from ((23l) . yl depends only on d and 7. The scal- 
ing function Hq depends only on 0, while the argument 
is scaled by the characteristic displacement v* multiplied 
by the scaled temperature Tm^ . This indicates that the 
scaling function of the thermal boundary layer exhibits 
universality for > since only scales of order Xjm con- 
tribute to the final scaling function, all features of g(u) 
with u ^ u„i being subdominant. 



5. Non-analytic cusp from the T ^ limit of the TBL 



The TBL functional (|134p exhibits a non-trivial large- 
argument limit, V = v/T —>■ cx) (or, equivalently, Qc — > 
00): 



ri[v]Ki—^ dqu{q)^/q~~q, (146) 



using that {\z\)z — •\/2/7r. This limit must match the 
T = limit of the functional R[v], denoted Rt=o[v]. 
More precisely, as we will show in the next subsection, 
ri must match the cubic term ~ v'^ of Rt=o[v], i-e., the 
cusp non-analyticity. 

We now perforin an explicit calculation for a uniform 
Vx = V. From (I142|) one finds in the limit of large vL'^^^: 



L'^TZiiv) = 7^lH w p(mVLV4)''/', 



Gc 



dGu{G)y/G{uc)-G, 



(147) 
(148) 



G„ 



which produces a cusp non-analyticity at zero temper- 
ature. The coefficient p can be computed explicitly for 
models I and II using Eq. ()143|144p and one finds: 



16 



— 9 — f) —9 — \ -^/^ 



(149) 



X2F1 



3 

2' 2" 



0'2' 



1 



/■mc\ 
\ m ) 



2+0 



One can check that in the limit m <C vric^ one has 
p ^ m^('^+^/2) and that the prefactor agrees, using (|147p . 
with the one obtained from the large v limit of (|145p . We 
have thus obtained the exact leading non-analyticity at 
zero temperature: it consists in a linear cusp in the force 
correlator, —N~^^j^d'^.R[v\ ~ |f;|. This non-analyticity 
is of the same kind (i.e., proportional to \v\) as the one 
found in the standard FRG (e.g., to one loop for any iV), 
and in the large N FRG in the regime v^ ~ iV, which 
will be discussed in detail in the next Section. However, 
there are differences in the dependence of the amplitude 
of the cusp on 771, L and A'^. Anticipating the results 
of Sect. IIVI they can be summarized as follows. In the 
regime m <^ rric one finds: 



RT=a[v] 


— Rt= 


oM^L^ 


Rt- 


-o{v) - 


- Rr^oiv) 


(1 


Gi 

m29 


r vV 

yvj 


'l-^M + o 




■■--AS 


Cn 


r vV 


1 — a 

Vc 


f 




, v^=0 


iNL° 
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where a,b, as well as Ci,Cjv are numerical prefactors 
(with Ci/Cn = Gc/G,n > !)■ Note that the two 
scales w* == {mL)^'^i'^mr'^ and Vc — \Nm~'^ are very 
different. The non-analytic corrections in the regime 
v^ — 0{L~'^/'^N'^) are related to the occurrence of shocks 
in the system, as the quadratic well is moved. Their con- 
tribution starts dominating once w > w*. 

From ()149|) one finds that the amplitude of the cusp 
is zero at the Larkin mass nic, and then grows linearly 
as function of nic — rn. This is in contrast to the large-w 
regime, where the amplitude jumps to a finite value at 
rric, cf. Eq. ((^M)) below. 

Let us finally point out that a cusp non-analyticity 
proportional to_|w| in the regime w^ = 0{L~'^/'^N^) was 
found in Ref. |47J , but with an amplitude scaling differ- 
ently with 771. This resulted from a calculation of a rather 
different observable which is reviewed in Appendix [E] 
A closely related cusp singularity was also found in the 
study of shocks in Burger's turbulence. ^^ (for a recent 
discussion of their relation to FRG see Ref.'*^). 

It is important to note that the non-analyticity found 
here at T = in the regime v — 0{\) is a robust feature 
that occurs irrespectively of the type of the ultrametric 
RSB scheme (whether continuous or 1-step, marginally 
or fully stable). As will be discussed below it reflects the 
switches in the minimum energy, and shock- like jumps in 
position, which occur as the energies of two states cross 
upon moving the harmonic well. It is, not surprisingly, 
rounded by temperature. In that sense it has some sim- 
ilarities with shocks discussed for interfaces A'^ = 1, and 
for Burgers turbulence. In contrast, in Section IIVI the 
cusp in the regime v^ — 0{N) will be seen to rely on 
the marginality of the RSB-scheme with respect to clus- 
tering fluctuations on the largest scales. That type of 
marginality only occurs naturally for systems exhibiting 
continuous RSB. As we will see, in that case a cusp oc- 
curs even at finite T, in contrast to the above discussed 
non-analyticity which forms only in the limit T = 0. As 
we will discuss later, it is related to a more global trans- 
formation of the set of states as v is varied. 



6. T ^0 limit ofR[v] 

The above described perturbation expansion is also 
perfectly suited to analyze the limit T — 0, where it turns 
into a rigorous expansion in |w|. A similar expansion was 
pointed out in Ref. |43 . To exhibit the dependence on 
T and v, we define q — v'^L'^j/T'^ by introducing the 
reduced coupling function: 



7(u) 



1 



9k' 



Vk 



I Id 



G(fc,u =Tu) 



(151) 



where in this paragraph we denote by v^ a suitable av- 
erage of Vk, such as v'^L'^ — J^v^ = Jk^k^ ^^^^ that 
Vk/\v\L'''^^ is just a form factor. As mentioned before, 
7(u) and its inverse 0(7) have finite T = limits, taking 



values between j{uc,m) — lc,m and Ucm, respectively. 
After rescaling y = \v\L'^/'^y/T^ the evolution equation 
([M]! for M{q,y) = M{'-f,y) becomes: 



d,M = -- 


'dlM+\v\L^/^ui^)dyiM^) 


, (152) 


M(u„y) - tanh (£!!^\ ^^0 sign(y), (153) 


md the second cumulant takes the form; 


R[v]^R[0] ^=" -\J gk^GT=oik,Um)\vkf (154) 


+2 


v'^L'^ dyy[M{u,rr,y)-sign{y)]. 



From (|152p it is easy to see that the above described 
procedure corresponds to a power series expansion in \v\, 
the term corresponding to 7?.„ "^ [v] being proportional to 
|t;p+". We will show below that in the case of a one step 
RSB the zero temperature correlator can be obtained in 
closed form. 



7. Comparison with the TBL from droplet arguments 

It is interesting to note certain analogies with formu- 
lae obtained from droplet arguments, where one assumes 
rare quasi-degeneracies of the ground state. In d = it 
was found in Refs. lisHiol that: 



B':Av):^dl,,,Riv)=R-M + 



(155) 



m T{y'y^F (m y — 



o(r'), 



F{z) =. ^ coth(z/2) - 1/2 = i [ij (z/2) - 2] ,(156) 

where y = u^^ = u'^ — u^ is the difference between the 
center of mass displacement of the ground state (u^) 
and of the excitation (u^). They are characterized by 
an (unnormalized) distribution D(y) = P{y,E = 0) 
where P{y, E)dydE is the (normalized) joint distribution 
of position and energy differences between the two states. 
Here we denote (...)j, := J dy...D{y), which can be nor- 
malized using the STS identity {yiyj)y = 2dij/m'^. In 
some simple cases (e.g., the Sinai model, corresponding 
to A'^ = 1 with random field disorder) D{y) is known an- 
alytically. The above formula can be generalized^!^ to 
any d: 

with J^" {z) = F{z), and a formula for all cumulants was 
also obtained. 
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Given the similarities between the formulae (|155p and 
(|134p . it is tempting to interpret the latter in terms of 
a droplet size density, summed over all ultrametric dis- 
tances u. To this end, we rewrite (|134p for uniform 
V = v/T as: 



2rp3 



L^Uiiv) = rr 



dqu{q){^P{\v\^2{q,- q)z)), il58) 



where we note that g(u) := q{u)/v'^ and its inverse func- 
tion are independent of v. 

As shown in App. |F1 the force correlator splits into a 
longitudinal and transverse part, describing forces par- 
allel and orthogonal to the displacement vector v. Note 
that due to the 0{N) symmetry the correlator is only 
a function of \v\. The TBL contribution to the force 

correlators is L'^A^^^v) = -L''7^i'(|^;|) and L'^A^r^\v) = 
—L'^Tl[{\v\)/\v\, respectively. They can be cast in the 
form (for details see App. IFJ) : 

/>oo 

-L^A^l\{v)=T dbp^f{b)^P{b\v\), (159) 

with the densities 

pT^'ib) = (160) 

"' dqu{q) ((z^ - z^)S{b - zv/2(gc - g)))^ , 
p?SB(^) ^ (161) 

''' dquiq) ((1 + z^)S{b - z^2{q, - q))) ^ . 

This has precisely the same form as the droplet expres- 
sions for the force correlators, which from (|155p are found 
as: 



» drop 
~'-^L,T 



iv) = 



(162) 



i?"(0) - m^T + T dbpL°r^{b) ^{b\v\) + 0{T^ 



where i?"(0) = d^. R{v = 0) (for any fixed i), and the den- 
sities P]^'^ (6) are given in (jFSp . Note also that formulae 
P32ll33p . and the appearance of the function f/;, bear 
similarities to expressions obtained in the case iV = 1 
for averages over a uniform density of shocks rounded by 
temperature^!^. 

In order to go further in the comparison and extract 
observables such as shock and droplet density and their 
size distribution, a more detailed description of higher 
moments is needed for the present model. Work in this 
direction is in progress. 



[= Ui in the notation of ([89)1 ] and the two self-energy 
parameters ctq,! (I26p . Similarly, the correlation func- 
tion G(fc, u) and the coupling q{u) assume only two off- 
diagonal values Go,! and qo.i- The former has the follow- 
ing interpretation: The measure on configuration space 
describing fluctuations both due to disorder and thermal 
noise can be constructed independently for each mode k, 
following Ref. \M- 

For the displacement in each environment one picks a 
set of "state centers" w^ as2^ 



uI + ^gM^ 



G,{k)^^, 



(163) 
(164) 



where the ^^ , ^^ are chosen from independent univariate 
Gaussian distributions. Note that in a given environ- 
ment the (infinite) set of m^ are globally displaced by 
the same (random) u° . Each state corresponds to a par- 
tial Gibbs measure which is assigned a weight Wa (with 
J2a ^a = l)i such that the total Gibbs measure in this 
environment is the weighted superposition of the partial 
Gibbs measures, i.e., in one thermal realization one picks 
a state a with probability Wa and the mode Uk takes the 
value 



Uk 



+ ^G{k)-G,{k)vt , (165) 



where 77^ are chosen from independent univariate Gaus- 
sian distributions and account for thermal noise inside 
a given state. In each environment the weights Wa are 
independent random variables chosen from the distribu- 
tion P{W) '--^vy^oo M^~(^+"'^\ the glass transition corre- 
sponding to the divergence of its first moment at large 
W , which implies that only a few states dominate the to- 
tal Gibbs measurciS^ An analogous construction applies 
in the case of continuous RSB, where the obvious gener- 
alization of the above to K-siep RSB has to be taken to 
the limit X ^ 00. 

With a general 1-step Ansatz of the form P^ we ob- 
tain the correlation functions: 



Go(fc) = Taogl, 



T 



Gi(fc) = TCTofffc + — I gk _i , ^ 



1 



(166) 
(167) 



where Si :— [(t](uc) = Uc(cri — ctq) As always, on the 
saddle point, STS implies G'^{k) = Tg^, and we note: 



D. The case of 1-step RSB (d < 2, 7 > 7^) 

As mentioned before, the GVM saddle point for a 1- 
step situation is characterized by the "break point" Uc 



G{k) = G%k) + UcGo(fc) + (1 - Uc)Gi(fc). (168) 



With such a 1-step Ansatz the GVM free-energy density 



18 



takes the forni)^ 



1 

+ 7: 



T. 



5,T^G(fc)--TrlogG(fc) 



(169) 



fo + ^[ucB{xo) + {l-Uc)Bixi)] 



(170) 



+ - 



Tl-Ur 



Si 



-In 



fc^ + m^ + El 



where we denote the square of transverse intra/inter- 
state fluctuations as: 



Xo,i 



G(fc)-Go,i(fc) 



(171) 



and we have absorbed quantities that depend only on 
T and m^ but not on Uc and Si, into the constants /o 
and /q. The free-energy density p69p is minimized with 
respect to G by the saddle point satisfying 5T /5Gab = 
[Eq. p^ ]. taking the 1-step form: 



CCi 



T 



B'{xo,i)- 



(172) 



In order to describe equilibrium thermodynamics, the 
breakpoint Uc has to be chosen so as to extremize {max- 
imize) T , as usual in replica treatments. However, other 
choices of Uc are of physical interest as well, as discussed 
in detail below. 



1. Instability of the RS solution 

Let us derive the phase diagram in the 1-step case. For 
simplicity we restrict to d < 2 and K/rric = oo. We recall 
that we use the natural units introduced in Section [TTl 
Performing the integrals, the free-energy density reads 
(dropping the constant /): 



H^c) 



1 



NL'i 2T 
AdT 1-u 



u,S(xo) + (l-u,)S(xi)] 
Si 




(173) 



where we have used 
1 



1 



r(i - d/2) 

(47r)'^/2 



Ei)2 -m" 



2A, 



e(2 - d) 



(174) 



In d = (a particle) the last line of (|173p becomes 
— 21n(l-|-Si/m^). In 0(uc) above. Si and Uc can be con- 
sidered as two independent variational parameters. The 



variation with respect to Si (at fixed Uc) yields back the 
saddle-point equation: 



Si = -^[B'(xi)-S'(xo)], 



(175) 



Xo 



Xi 



1 (|17ip one has: 
4Ad T I 1 


1-Ue 1 


e(2 - d) Uc 
4Ad 


T 


(m2 + Si)i-'i/2j 



e(2-d)(m2 + Ei)i-'i/2' 



(176) 



The replica-symmetric solution (uc = 1, Xi = Xo = X^^i 
Si — 0) is vahd at high temperature/large mass, but 
becomes unstable when the condition 






1 



(177) 



is met, with x^^ = £(2-'^') -^"^'^ ^- This defines a unique 
function Tc{m) for masses m < rric smaller than the zero 
temperature critical mass:^^ 



rric = mc{T = 0) 



4^„ 



B"{0) 



i/e 



(178) 



The function Tdjn) describes the location of a contin- 
uous glass transition towards a 1-step RSB phase for 
rric > JTi > m^. Here, to* denotes the mass where 
Tc{m) attains its maximum. For to < m^, the line 
Tc{m) has little physical significance, since in that regime 
the glass transition occurs in a discontinuous manner at 
T > Tc{m), as will be discussed below. 

It is interesting to note^^ from (|163m65p . using (|166I 
11681) , that as TO ^- the thermal fluctuations of the dis- 



placement field u within a state, G ~Gi — T jik^ + to^ + 
Si) remain bounded and massive and occur only at the 
Larkin scale. 

By contrast, the r.m.s. difference in displacement be- 
tween two states in a given sample, m" — w" , scales as 
Gi — Gn = — [ i.a i"^ a — -2 I "^:j I ^ ] . Hence, the ther- 
mal fluctuations between states, which are active at any 
T > in a given sample, occur at all scales up to 1/to. 
Note that these r.m.s. differences in displacement have 
a finite T = 0^ limit, if T/uc goes to a constant in that 
limit (which is generally the case as shown below). These 
fluctuations are responsible for the roughness of the man- 
ifold C = (2 — (i)/2 with a non-zero amplitude even as 
T^O. 

Finally, sample-to-sample fluctuations include in addi- 
tion the term Gq ~ rao/(fc^ -t- m^)^ whose mass depen- 
dence is controlled by ctq- These fluctuations also occur 
at all scales smaller than Xjm. 



2. Metastable states and configurational entropy 

The break point < Uc < 1 is a priori a free parameter 
of the 1-step Ansatz. In order to describe the thermo- 
dynamic equilibrium, one should choose Uc = u°'^ which 
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extremizes the free energy: 
0, 



^K.') 



(Equilibrium) 



(179) 



(at fixed Si) as was done in Ref. |4y. As in most glassy 
problems with a 1-step structure, at low temperatures 
T -^ Tc the equilibrium break point is proportional to 
the temperature^* uj;*^ « Tlu°'^. This is shown in App. ]G\ 
where a;°'^ is computed as well. 

However, the function 4>{Uc) encodes much more in- 
formation than just the equilibrium physics. It can be 
interpreted^ as the quenched "replicated free energy", 
i.e., the free-energy density per replica of a set of Uc 
clones (replica) which are bound to remain in the same 
metastable state of the energy landscape. With this in- 
terpretation one can alternatively write the total free en- 
ergy as a sum over states. Their number Af{f) at fixed 
free-energy density / increases exponentially with the 
volume and N. One thus introduces the configurational 
entropy, or complexity, E(/) = log[U {f)]/NL'^. The to- 
tal Boltzmann weight of Uc clones can then be written 
as: 

exp [-NL''puc(biuc)] (180) 

«yd/exp{-7VL'^[/3u,/ -£(/)]} , 

where /3 = 1/T. In the large- A^ limit, the configurational 
entropy becomes the Legendre transform of the replicated 
free energy: 

/3u,<^(u,) = /3u,/ - S(/), S'(/)=/3u,. (181) 

Knowing '/>(uc), one easily obtains an implicit 
parametrization of the configurational entropy as a 
function of the free-energy density: 



d[uc 



duc 

S = /3u20'(u,). 



Uc0'(Uc 



(182) 
(183) 



We see in particular that the equilibrium prescription 
(|179p corresponds to choosing the states with vanish- 
ing configurational entropy. Since the configurational 
entropy is an increasing function of /, the so selected 
states have the lowest free-energy density in typical sam- 
ples and thus describe indeed the quenched equilibrium 
free energy. 

However, the equilibrium may not be possible to reach 
dynamically, which suggests a different choice for Uc- The 
choice is constrained by the requirement that the 1-step 
solution be stable (as was the case for Uc = u°'^). 

Metastable states of the free-energy density above the 
equilibrium value f°'^ are described by breakpoints in the 
range u*'^ < Uc < u'^'^, where the lower boundary is deter- 
mined by the so-called replicon instability of the 1-step 
scheme, i.e., the conditioni^S 



4A, 



A comparison with (|177p shows that at fixed temper- 
ature T < Tc{m^), this impUes m^ -f Ei = m'^{T), 
where mc{T) > m* is uniquely defined as the solution 
of Tc{mc{T)) — T. One can prove that solutions of the 
1-step saddle-point equations with condition (|184p exist 
for all T < Tc{m.^) and m < mc{T) (but nowhere outside 
this parameter regime). ^^ The states described by u*^ 
are merely marginally stable and are often referred to as 
"threshold states" . Since they are usually the most abun- 
dant metastable states of the system they are most likely 
to trap the dynamics after a fast temperature quench for 
m < m^,. 

Another choice for Uc of interest is the value Uc — 
u^P > u^'^ where the 1-step scheme becomes unstable 
with respect to a clustering fluctuation among the ex- 
isting states^^. The latter is equivalent to the condition: 



4A„ 



B"{xo) = 1- (Cusp for FRG at large v) (185) 



This condition ensures that the second cumulant in the 
regime v'^ ^ N is non-analytic at w = 0, as we will discuss 
in Sec. IIVB4I 

We emphasize that the kind of marginal stability im- 
posed by (jl85p is clearly distinct from the marginality 
P84p due to the replicon mode, which is usually imposed 
to select dynamical threshold states in 1-step systems!^. 
The latter is also the marginality property of continuous 
RSB systems that ensures the presence of collective soft 
modes in classical^i. and quantum mean-field spin, elas- 
tic or electron glasses^^^i^i^, and leads to a universal, 
saturated pseudogap in the local field distribution of spin 
and electron glasse o^^'^^i^^ . 

While the replicon instability (I184p indicates the 
fragility towards fragmentation of existing states into 
substates, the condition (|185p signals the instability to- 
wards the formation of clusters among previously equiva- 
lent states. When m <C rric these two instabilities involve 
rather different length scales. This can be seen as follows: 
The replicon instability corresponds to each state u^ in 
()163p giving birth to a new cluster of substates labelled 
by /?, u^ -^ u^ = u^ -)-(Su^, where each Su^. occurs with 
a probability W'a (with ^o W'a — 1). This rearrangment 
implied by the additional substructure involves scales of 
order 1/toc since dG{u) — Tda{u)/{g'i^^ + [o-](u))^ with 



[a](u > Uc) = Si 



On the contrary (|185p signals 



the instability towards u^ splitting into several clusters 



e(m2-hSi)2 



B"{xi) = 1. (Threshold) (184) 



u^ ^ u^ + 5u^, generating substates u^'"- This rear- 
ranges the original u^ into clusters labelled by /?, and 
obviously involves scales up to 1/m. 

The above two kinds of instabilities have been analyzed 
in detail in the context of the spherical p-spin model 
in Ref. iZS. More recently they have been discussed in 
the context of lattice-glass models^ and optimization 
problems^SiliiiZ^^ where they are sometimes referred to 
as instabilities of the first (clustering) and second (frag- 
mentation) kind^"^. They also appear as the instabilities 
driving the transitions between the two types of one-step 
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FIG. 4: Phase diagram for model I in d = 1. The solid line 
for m > m* indicates a continuous glass transition, where the 
renormalized force correlator B'{x) displays a cusp at the ori- 
gin. For smaller mass, m < m*, the glass transition towards 
the one-step RSB phase (as a function of T) is discontinuous 
and takes place dynamically as a freezing transition at Td, or, 
if equilibrium can be attained, as a genuine thermodynamic 
transition at the lower temperature Tc- A similar phase dia- 
gram applies in d < 2 for model I and model II with 7 > 7c. 
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FIG. 5: The breakpoint Uc as a function of m at constant 
temperature T = 0.8 in the glass phase (m < rndT) = 0.673) 
for model I in d = 1. The labels cp, eq, and th indicate thresh- 
old, equilibrium and cusp states. Notice that u'^^ exceeds 1 
at small enough m. The same would be true for u'^'' and u**^ 
for r™'"' > T > Tc(0) and T™"" > T > Td(0), respectively. 



on the instability line Tc{m) where (I177P is obeyed, there 
is only one admissible value for u^: 



phases adjacent to an intermediate two-step RSB phase 
in certain mixed spherical spin modelsi^ 

For model I in d = 1, i.e., the directed polymer prob- 
lem, we have checked explicitly that u**^ < u^'^ < u^p, and 
we expect this to be generally true for 1-step solutions. 
This is illustrated in Fig. [S] For the associated config- 
urational entropies this implies E**^ > S'^'^ = > S'^p, 
and thus the states selected by the clustering instability 
()185p have negative configurational entropy. Anticipat- 
ing the analysis of Sec. IIV B 41 we conclude that samples 
exhibiting non-analytic shock-like behavior in the regime 
v^ ^ N where an FRG equation was previously derived, 
correspond to exponentially rare realizations of the disor- 
der occurring with probability P ^ exp [— iVi'^|E(u^P)|] , 
as derived explicitly in App. [G) 

Note that upon approaching the limit of a marginal 
one-step solution, i.e., for d < 2 and 7 i 7c, the three 
values u*^ < u^'' < u^p merge and the one-step saddle 
point becomes simultaneously marginal with respect to 
both instabilities discussed above. 



3. Phase diagram 

In this section we establish the phase diagram and dis- 
cuss some of the subtleties associated with the choice of 
Uc. The situation is closely analogous to a particle in a 
random environment (d = 0) analyzed in Ref. UM The 
phase diagram of a typical case, model I in d = 1, is 
shown in Fig. [4l and explicit values given in the analysis 
below refer to this specific case. 

The glass phase exhibits everywhere 1-step RSB. Right 



Hm) 



4A, Tc(m)i3'"(x«S) 

£2 ^2-d B"{X^S) ' 



(186) 



which yields the fluctuation-dissipation ratio relating re- 
sponse and correlations at the glass transition. ^^ It in- 
creases monotonously with decreasing m, and reaches 
uj;"* = 1 at TO = m* (= 1/e = 0.3678, with T^^"" = 
Tc{m^,) = 3/e = 1.1036 in d = 1 for model I where 
Ai = 3/4). For states to be dynamically or thermody- 
namically of significance, u must always be smaller than 
1, and thus the line Tc{m) looses its significance below 

TO*. 

For m > TO* the glass transition is continuous in the 
sense that Gi — Go smoothly goes to as T t Tc{m). For 
TO, < TO* the temperature-driven transition is discontinu- 
ous with a sudden jump between intervalley and intraval- 
ley correlations occurring at some T G [Tc{m),T"'^'^^]. 
The location of the transition depends on the point 
of view. From a thermodynamic (static) standpoint 
the glass transition takes place at the line given by 
ul'^{m, T) = 1, indicating that the replica symmetry must 
be broken spontaneously to extremize the free energy, 
which gives rise to a non-analyticity in the free energy. 
For model I, one finds the explicit result 



Tc(to = 0) 



2Ad 
de 



{2ef^'' 



' d=\ 



0.9028. (187) 



However, metastable states exist already at higher tem- 
perature. In the mean-field limit A^ — > 00 they induce 
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a dynamical freezing transition at the line Td{m) defined 
by u*''(to, Trf) — 1. For model I one finds 



Tdim - 0) 



e2^/^Ad^ ''"'- 



d=\ 



1.0268. (188) 



The resulting phase diagram is naturally very similar to 
the one for a particle in a random environment (the limit 
d = of a random manifold)^!, and it also strongly re- 
sembles the one of the spherical p-spin model, whereby 
the mass m takes the role of the external field />2^. In the 
latter the glass phase is known to be everywhere of 1-step 
nature, the transition (at fixed K) being continuous for 
\h\ > h.if and discontinuous for \h\ < ft,*. 



J^. Second cumulant 

To evaluate the second cumulant of the renormalized 
disorder potential V, we introduce as before the coupling 
go,i, connected to Go,i via: 



90,1 = -^ / 5fc^Go,i(fc)|wfe| 



(189) 



We need to apply the recursion relations (|M|) only once 
(since here K — 1) to calculate: 



9i{y) = <2xp[uctlJi{y)] 



(190) 



= e^^^l^u. ( [2 cosh(2/ + zVgi - go)] ""'), , 



and using similar steps as in the derivation of 
(I103P we obtain: 

R[v] - R[0] = -2T'^qo 



and 



(191) 



~2T" 



dy {"01 (y) -ln[2cosh(y)]} 



Recalling that q'^ = q — [Uc^o + (1 — ^c)qi\ = 0, one can 
see that V'i(y) — ln[2cosh(2/)] ~ q'^ /2 = at large \y\, 
and thus the integral in (|19ip indeed converges. In the 
following we examine various limits of this formula. 



5. Large L''v^ 

Let us introduce the notation Q = qi — q^. At large v 
(large Q) one can use that ((2cosh(y + z^/Q))"'')^ ^ 
2cosh(uc2/)e^<^'5/^, as can be shown using a saddle-point 
method. This yields the large-u limit: 



R[v] - R[Q\ " -°° -2T^qo 



(192) 



1 f°^ 
~2T'^— / d2/{ln[2cosh(yUc)] -Ucln[2cosh(y)]} . 

Evaluating the integral we find the final result: 



T f TT^T^ / 1 



i?[0]-i?H"^°°-a(0) n 



6 V u 



-^ - 1 , (193) 



which matches the full RSB result (|107p . if we formally 
replace Um by Uc- These two expressions are indeed iden- 
tical in the case of a marginal 1-step solution, as it oc- 
curs in d = 2. In that case, 9 = 0, and the analysis 
of the crossover to large v in Section pil C 1[) remains 
unchanged. 



6. Small L''v^ 



For small L v one can expand in Q as follows: 

V'i(z/)-ln[2cosh(y)] (194) 






(l-u,)[l-tanh2(y)]Q + 



where q'^ = was used. Performing the j/-integrals, one 
eventually finds: 



R[v] - R[0] = -2T^qo - 2T2(1 - uj x 



(195) 



Q-^Q' + ^Q' + |^(3^c - 5)Q4 + OiQ') 
6 45 315 



As it must be, the leading term is: 



R[v] - R[0] = -2T^ / qiu)du 



(196) 



which yields the GVM result for two-point correlations. 
Indeed, (|196p matches the expression (|113p . which should 
not depend on the RSB scheme. Similarly, inserting the 
one-step form for q{u) in (|114p one checks that it repro- 
duces the second-order term from (|195p . In fact, one can 
check that (|195I114|) agree with: 



R[v] - R[0] = ~2T^ 



tr(g) + -tr(q2) 



(197) 



-±r£ql^ + 2tiiq')-3ty{q'Mq)\ 



315 



(-3tr(g4) - 12[tr(g2)]2 + 30tr((j2)[tr(<j)]^ 



+5^g^^-20tr(g)^g3^ + O(g^) 



where we have used q'^ = and defined tr(^) = 
lim„^o ^Tr(A). Note that J2b1ab = HiQ ' q)q), where 
the dot is the Hadamard-product. 



7. Low-T expansion 

Let us now consider the low T expansion, i.e., the 
thermal boundary layer, using similar notations as in 
PIIC3|) . In the T ^ limit u^ -^ 0, but Uc = Uc/T -^ 
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u)'^'^ has a non-zero limit, and so do go and 51 when ex- 
pressed in terms of fixed ii = v/T, as in (|116p . Hence 
it is convenient to define, for the present 1-step case, an 
expansion in Uc^-^, analogous to (|12ip : 



R[v] 



p=0 



(198) 



obtained by expanding (|194p in powers of Uc at fixed Q. 
This yields: 



7^oM --2rV, 

/oo 
dy 

— /in cosh (?y 



lncosh(y + z\jQ) 



Q 



Q) 



(199) 



(200) 



(201) 



dy ( lncosh(j/ + z-\/(5) — Incosh(y) 



p+i 



The expression for TZq has been simplified using an inte- 
gration by parts and the identity ()102|) . and is found to 
match the result (I126|) for the continuous case. The term 
q — qi ^ ~UcQ has been included in TZi. The term mul- 
tiplied by u^ in (|20ip is the {p + 1)**^ connected average 
(cumulant), as indicated by the superscript c. After a 
calculation summarized in App. iHl we finally obtain the 
thermal boundary layer: 



rQ I 
7^lH = 2T^Uc / dq{z^qcoX\\.{z^/2q) 



(202) 



which agrees with the result for continuous RSB (|134|) 
as expected, since both formulae should apply for the 
borderline case of a marginal 1-step solution. 



8. Cusp and full correlator atT = 

The case of one-step RSB is sufficiently simple to allow 
for a complete calculation of the non-analytic T — limit 
of the full functional R[v] in the regime v"^ ^ 1. Let us 
define the variable: 



w := 



Q 



gf[G^ik)~Go{k)]\vk\ 



1/2 
(2Q3) 



which will be shown to be of order 0(1) when the 

crossover to the shock-dominated regime occurs, see 

((2T2| below. It tends to a finite hmit as T -> 0. For 
a uniform v^ — v this becomes: 



^^-[L'uym\G^~Go)f' 



V 

2^ 



Si 



Si 



1/2 



(204) 



with Go,i = G'o,i(fc = 0) as given by (|166ll67p . Note that 
the characteristic scale for v (determined by w « 1) is 
not exactly the scale w* — L^'^^^m^^^^ for shocks found 
for continuous RSB in p09p . but rather (for m <C Si) 
~ w*** "^^ = v^,m2 r^ L^'^i'^m^^ . However, the two scales 
become the same in the case of marginal one step RSB 
since there 9 = Q. 

Let us first obtain the T — Q cusp for a uniform Vx — v 
by taking the limit of large v/T of the thermal boundary 
layer. In that limit the expression (|202p tends to^S^ 



L'i7^f=0(«) = 2u 
3(7ru,)i/' 



itf 



Si 



Si 



3/2 



1 8u;3 
(205) 



using {\z\)z = \flpK. This finite non-zero limit as T — > 
(cf., App. lGj) is thus the exact expression for the leading 
non-analyticity of R\o\ — _R[0] at T = 0, which is again 
cubic '--^ |vp (corresponding to a "linear cusp" of the force 
correlator —B!'\v\). Note that in this regime, v^ ~ i""*, 
the cusp exists irrespective of the choice of Uc, including 
the equilibrium solution. This is to be contrasted to the 
regime v'^ ^ N where, in the case of a non-marginal one- 
step solution, a cusp exists only for the choice Uc — u'^p, 
see the discussion in Sec. IIV B 41 

We now obtain the full functional Rt=o [v] for arbitrary 
Vx, using the variable vu defined in (|203p . By using (|190p 
in p9ip and taking the limit of small Uc, we find at T = 0: 



R[v]-m\T=o^-l 



g-'G'o^^>{k)\v,\ 
+TlT=a{w), 



(206) 



where: 
TZT=oiw) 



= -7^ / dy 



8ydy 



l+crf[(u.-j/)/V2] 



8H 
3v^ 



0(w^ 



w< 1, 
w> 1. 



(207) 



From the expansion for w ^ 1 we again recover (j205p 
with, in the uniform-u case TZt=o{w) -^ L'^Tif^'^{v). The 
asymptotics for tu 3> 1 is in agreement with the result 
(|193p for continuous RSB (with Um -^ Uc). 

From this expression one can compute the force corre- 
lator, which, from 0{N) symmetry, splits in transverse 
and longitudinal parts as defined in App. [Fj 

The non-analytic parts of the longitudinal and trans- 
verse force correlators in the regime v"^ ~ L^"^/"^ are 



proportional to 



-T^T-- 



q(w) and 



Ky^o(tu) 



, respectively^ 



They are plotted in Fig. [SI both exhibiting a linear cusp 
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FIG. 6: Nonanalytic contribution to the force correlator in 
longitudinal Al{w) and transversal At{w) direction, at T = 
in the case of 1-step RSB. The rescaled displacement w is 
defined in (l204ll. 



at the origin. Note that the full correlator of the force in 
addition contains a constant which derives from the v'^ 
term in ()206p . Further decay from this constant to zero 
at infinity occurs (for any 7 > 0) in the regime v'^ ^ N. 



E. Interpretation of u, and the thermal boundary 
layer 

What happens to the manifold as the center of the har- 
monic well is moved? Let us first focus on the case of a 
one step RSB for which the construction of the statisti- 
cal measure is significantly simpler and was recalled in 
P63ll64p . From this construction we know that at low 
temperatures a certain state center S.^ dominates, but 
there is at least one excited state which differs in energy 
by a typical amount T/u^,. Changing the position of the 
harmonic well by an amount v shifts the relative energies 
of these two low lying configurations (states), noted u^'^, 
by an amount: 

AE^ I g'^^Vk-[u\k)~u^{k)]. (208) 

Jk 

One can expect a switch from one state to the other, and 
thus a macroscopic jump, as the energies of two states 
cross. This occurs typically when \AE\ ^ T/uc — u^^. 
Recalling that [u^{k) — u'^{k)]^ is of the order of Gi(fc) — 
Go (A:), one finds: 



AE^^2 g^'\vk\'[Giik) - Goik)]. 



(209) 



Hence, comparing with Eq. (j203p the jump occurs typi- 
cally when: 



u;-0(l). 



(210) 



We thus understand the scale of variation of the func- 
tional R[u] given by (|206l207p . The argument is even 



simpler for a uniform center v {vk = L'^^'^v6k,o) which 
couples only to the A: = mode of the excitation, 
the energy shift being to^L'*/^u[u^(0) — v?{Q)\. Since 
[ui(0)-u2(o)]2 ^ 2(Gi - Go), we see that the energy 
shift (|208p competes with the typical excitation energy 
u^^ when v is of the order of w*^ "'^ as defined below Eq. 
(|204p . Since the energy difference is finite, the jump is 
thermally smoothed at finite temperature, but turns into 
a sharp shock for T -^ 0, which is at the root of the 
non-analyticity in (|207p . 

We can obtain better insight into the physics for shifts 
of order v ~ w*^ °^ by examining the force-correlator. 
We consider uniform v for simplicity. To be specific we 
define —R"{v) := —d?.R{v) which corresponds to the 
longitudinal force correlator. From (|206p . we see that for 
V S> w*^ °^ the correlator is simply constant and equal to: 



R"{v > vl'^'P) = m^Goik = 0). 



(211) 



This has a simple interpretation in terms of the hier- 
archical construction (|163ll64p : One can imagine the 
states (characterized by ^") as parabolic potential val- 
leys within a big parabolic potential valley as determined 
by the global center u^, and accordingly write schemati- 
cally V{u) = VQ{u°) + Va{£,°') for the total potential. For 
V ^ v^,^ °^ , the manifold has jumped to a different state 
center, and therefore the part of the force arising from 
^q(C") does not contribute to the correlator. However, 
u° is very robust under shifts of v <^ N^'"^, and jumps 
only once v'^ ~ iV, as we will see in the next section. 
The force corresponding to the big valley, V^'(u°), thus 
remains constant. Since the disorder potential competes 
against the quadratic well L'^m?{vPY /2^ that part of the 
force is of order Vq{u) ^ L'^^'^m^u^^Q, which leads to 

dnn). 

For V <^ w*''*°^ we see from (I205P that the force corre- 
lator behaves as 



-R"{v<t:vl'''"P) 



m 



m^Gi{k^O) 



Istcp , ^ \l/2 



Si 
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3/2 



The first term can be understood as to* times the intra- 
state correlator (u^) . The force correlations rapidly 
decrease with growing v. Indeed, we can understand the 
non-analytic piece ~ \v\ as being due to shocks which oc- 
cur with a finite density along the w-axis. Note the large 
prefactor r^ L'^/'^ of |w|: It indicates that the product of 
the density and size of shocks scales as ~ L'^/'^, that is, 
in such a way that the system size rather than \/m acts 
as a cutoff. Further work is in progress to clarify the 
properties and the distribution of these shocks. 

Let us close the discussion of shocks by sketching the 
analysis in the case of continuous RSB. The typical scale 



in that case, w* 



Istcp 
— V " 



-e/2 



can be retrieved by the 



same argument as above for the one-step RSB. However, 
now the energy difference between two states at ultramet- 
ric distance u is typically of order T/u. This competes 
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with the energy shift induced by the displacement of the 
well: 



AE^ 



2/ 9f\vk\'[G{u,,k)-G{u,k)] 

Ik 



ai+2/e 



(213) 



The largest response is due to transitions between states 
differing at the largest scales, i.e., u Ri Um ^ mP ■ The 
above energy comparison at this scale immediately yields 
Vsf. However, there is no sharp selection of shocks that 
involve the highest hierarchy scale. Smaller shocks con- 
tribute as well, even though their weight decreases with 
increasing u characterizing the shock. Note that from the 
construction of the measure of pure states it is clear that 
when a shock occurs at scale u = Ur (in a finite if -step 
RSB scheme), the whole sub-hierarchy of state-centers 
u^ with K > s > r changes. Note that in general shocks 
characterized by a larger u correspond to smaller spatial 
rearrangements. 



IV. REGIME v^ 



N 



We now compute the functional W[v] = W[j = 
(g^^v)/T] defined in Section FlI C[ and more explicitly its 
second cumulant part R[v], in the regime w^ ~ N. This 
is the standard regime for the large- A^ analysis, and in 
Refs. [lllilSl the associated saddle-point equation was ob- 
tained. There the focus was on the calculation of the ef- 
fective action T[u\, and in particular its two-replica com- 
ponent R[u], part of which will be of use here, too, to 
obtain W[v] and R[v]. We find, as announced in Section 
IlICi that the functionals R and R are identical. The cal- 
culations carried out here are, however, rather different 
in spirit from the one in Refs. [lllllSl . We obtain R[v] 
from imposing external fields v"" = ±v to two groups 
of replicas, whereas in Refs. IlllIlSl the case of u"'' ^ 
for all a ^ & was considered. It is reassuring that the 
two approaches yield consistent results. The advantage 
of the present approach is that it includes RSB effects, 
and hence it provides a complete derivation of the FRG 
equation on all length scales (below and above the Larkin 
scale). 



A. General saddle point 

To analyze the saddle point in the regime v'^ ^ N it is 
convenient to introduce the rescaled variables: 



We can now rewrite the definition (|59|) and decouple the 
partition sum with the help of two auxiliary fields as: 



^W[v] 



S = 



11 z. 

a 


■a ^ y 

J J, 


^ 1 V[u]V[x]V[a]e-^^ , 


^}]f9k' U'^-k-ul^il.-ut) (215) 


^l 


u{x.) + ly 

ab 


<" (xf - K ■ «') 


: 



where the replica matrix field Xx'' has been introduced 
through a purely imaginary Lagrange multiplier matrix 
a'^'^. A standard choice, as in previous sections, is g'jl^ = 
k^ + m? but our analysis is more general. We have also 
included the source term shifting the center of mass. U 
is a function of a n x n replica matri: 



104 -,--, 



u": 



ab 



.-M2 



ur) 



1 



3!r3 



^5(3)(_ 



(216) 



ibc 



containing the information about the bare disorder via 
its cumulants, e.g., its second cumulant being (as in (l7|): 



Ra{u) ^ NB{u^). 



(217) 



Note that the rescaling (|214p was performed in order 
to make explicitly appear a factor of N in front of the 
action 5, which allows for a saddle-point analysis. One 
now explicitly performs the functional integration over 
the field Ux to obtain: 



W\v\ 



= / V[x\V[cj]c 



-NS[X:<tM 



S[Xi cr, v] 



— Tr In 
2 



yCr'-^) 



(218) 
(219) 



•'^ L ab ^^ 

9T ^-^ 



Vx [(.9 - .9^3) 



-liab ~b 

\XX' ^X' ■! 



where the inversions and the trace are performed in both 
replica and spatial-coordinate space. We use the short- 
hand g~^ = [g~^)'^y = 9xx'^'^^^ which is diagonal in 
replica space. 

Eq. piSp has now the standard form for a saddle-point 
evaluation of the functional W[v] =: iViy[-ii] except that 
the saddle point is not, in general, uniform in space. At 
dominant order in 1/iV, we obtain 



W[v] = ^lnJ2{e 



-NS[xi,<yi,v] 



(220) 



u 



V 

N 



(214) 



We have allowed for different saddle-points ( "sp" ) to con- 
tribute. To alleviate notations, we did not add an index 
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indicating the different saddle-points to Xv and an . They 
depend on v^ and are the solutions of the saddle-point 
equations obtained, respectively, by setting to zero the 
functional derivatives (at fixed Vx)'. 



SS[xv,(^v,v] 



Serf 

6S[x,<J,v] 





= 0, 


(221) 


X=Xi,c=<^ii 








= . 


(222) 


X=Xv,cr=t^v 







Sxf 



Solving these equations, one can obtain R[v] from its 
definition ^5^. 

The explicit form of the saddle-point equations is: 

1 



iXv):' = {G,)t + ^{G, : ff"i : v): ■ (G, : g"' : «)^ 

(223) 






ab 



2^2 ' 
\Gv)xx + '^x'^x ■ 

d 



-2T- 



g^ab 
TG^' = g-' - a. 



■UiiXv) 



(224) 
(225) 



where Ga is a matrix with both replica indices and spa- 
tial coordinates; thus inversion is carried out for both. 
The notation for the A^-component vector (G : j)^ := 
X^c/ ^xy^y ^^ ^ shorthand for a matrix product, but 
consistent with our above conventions, we do not write 
the vector indices of j explicitly. We have defined the 
average displacement induced by the source: 



1 



ul^ul{v):^^{G,:9-':i)l. 



(226) 



It is a function of v and of the chosen saddle point "sp" 
and to simplify notations we drop in what follows the 
dependence on v except when an ambiguity arises. As 
detailed in Refs. 111)13 , u{v) arises in performing the 
Legendre transform to obtain T\u] from Ty[i'], as one may 
see from (gH) and dUD, i.e.: 



(227) 
(228) 



T[u]+W[v] ^-u:g ^ : v , 



u^Tg 



5W[v] 
Sv 



and using that for a given saddle point: 

Wsp[v] = ~S[xd,crv,v]. 



(229) 



Taking the total derivative w.r.t v, which is equivalent to 
differentiating only the explicit v dependence in ()221|) . 
one recovers (I226p . While in Refs. 11 , 13 we had chosen 
sources v°' — v^ ^ such that m° — ■u'' ^ and therefore 
assumed a unique saddle point, here we allow for sponta- 
neous RSB and hence for many saddle points. In general 
the set of saddle points will include the set of all permu- 
tations IT which leave v invarian 
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I.e. 



y-^ia) 



all a = I, ...,n. Hence the u'^{v) in (|226p is identical to 
a partial average {u'^)tt (corresponding to a single saddle 
point labelled by tt) while the thermodynamic average 



corresponds to the full average over all equivalent saddle 
points, (u^) = J2Tri'^x)-^- I^ the hmit v ^ the sad- 
dle point equations (|223p become identical to the saddle- 
point equations of the GVM, with Gy — *■ G and (t„ — > cr 
taking the values of the solution given by Mezard-Parisi 
(see also below). Performing the sum over equivalent sad- 
dle points sp in (|220p one recovers the results of Section 

mil 



B. Analysis for a uniform v 

According to the general strategy to compute R[v], 
as described in Section IllICi we now solve the saddle- 
point equations (|22ip and (I222p . specifying the source 
v^ = Vx{l, 1, • ■ ■ , —1, • • • , —1) with n/2 entries -1-1, and 
n/2 entries —1. For simplicity we restrict here to a uni- 
form Vx ~ V, for which R[{vx = v}] — L'^R{v). From the 
definition ([5^ and the replica-sum expansion ([79]) one 
has: 



L-'^W[v] 



m 



,„~,2 



1 n^ 



2T 2r2 2 

R{v) = Ni3{v^), 



B{'ii^) + 0{n'^), 



(230) 



up to a (u-independent) constant. The matrices Xv and 
a are now independent of x, and Gn is translationally 
invariant. We parametrize the n x n replica matrix Xv 
by: 



Xi 



Xl X2v 
X2v Xl 



(231) 



where we anticipate that the diagonal blocks will turn 
out to be independent of v. We use similar notations for 
a and G: 



0'2l} CTl 



G5(fc) = 



Gi(fc) G2v{k) 
G2v{k) Gi(fc) 



(232) 
Note that ai and Gi should not be confused with the 
variables used for the 1-step solution. The saddle-point 
equations (|223p - (|225p now become: 



TG-,\k) 



Xl ^u^I+ Gi(fc), 

Jk 

X2v = -U^^ + / G2v{k), 

-0'2D fffc 11 - (Tl 



^ab 



(233) 

(234) 
(235) 
(236) 



and must be solved for a given value of v. We have in 

troduced 

matrices 



troduced the notation J = l„/2 and 11 for the § x ^ 



J 



ab 



1 , 



-na6 ^ cab 



(237) 



26 



In general, we expect that Xi is an ultrametric ma- of the external field v: 

trix while X2I = X2v^°'^, and similarly for a^'' and CTjc- 

This Ansatz for the solution is motivated by the ex- , 2 

■^ ab 

pectation that states in different groups will be very T 

distant in phase space (at least as distant as the far- 
thest equilibrium states, as described by u — 0), and 
thus their mutual overlap will not depend on the spe- 
cific replica in either of the groups. Further, we have 
u^ = 'u(l, 1, . . . , —1, . . . , —1), and one can show that our 
Ansatz implies 



B'ixr + xT-^xT) 



-5,,Y.B'{xT + xT-'^xT) 



ul = vl, (238) 

which means that the average displacement is tied to the 
center of the harmonic well. This holds because (|226p 
implies TgkUk = [GJ(fc) - Glj^{k)]vk, where G\{k) = 
Tib Gi,ab{k) (the sum being restricted to the n/2 replica 
in the same group as a) and G2y{k) = J2b G2v,ab{k) (the 
sum being restricted to the n/2 replica in the group 
not containing a). Within our Ansatz, as n — > 0, 
G1{k) = Tgt and G^-(fc) vanishes (see ([M5)) below), 
which establishes (|238p . 

To evaluate (I236P , let us split the set of indices a into 



for a, b in the + group. An analogous equation holds for 
the — diagonal block. This is not surprising since the 
external field v does not separate the replica within a 
group. Since these saddle-point equations have a unique 
physical solution, we will identify ai = a henceforth, 
and also drop the subscripts 1 on Gi and xi- In Parisi's 
parametrization of ultrametric matrices, the saddle-point 
equation is conveniently rewritten as 



a(u) 



--B' ( 2 
T 



(t(u) du. 



;G(fc) - G(fc, u)] 



(244) 
(245) 



for the first n/2 replicas, i.e. those for The saddle-point equation for the off-diagonal part is: 



two groups, 

which v°' — V, and a_ for the remaining n/2 replicas, 

and consider for simplicity a bare model with gaussian 

disorder-i2&: 



^2l = -j;B'{x 



aa I bb 

X 



2x2') 



(246) 



U{Xv] 



\ '' D/',,aa j^ bb ab ■.,ba\ 

^2T2 ^ ^-^1 -^1 ^-^1 ^Xi ) 

\ ^ aUA"- I ^,bb ^,ab ^,ba\ 

~2T2 ^ ^--^i -^1 ^-^1 ^Xi ) 



which within our Ansatz ctj^ — o-2vI"''' reduces to a single 
equation 



i_,b_ 



a2,=-^B'{4v^ 



;^ E ^(^r + Xi'-X2'-X2".)- (239) 



2 /[G(fc)-G2i(fc)]) . (247) 



2^2 



a_ ,64 



One can easily see that Eqs. (|236l239p imply the STS 
property cr"^ = Ta'^ii'' ~ ^- ^PPlying the inversion for- 
mulae p4j|36p to (|235p . we find the exact relation (to all 
orders of n) G^ — G^ + G2J5 = Tg^ from the connected 
part (|235p . Further, to lowest order in n one finds [cf., 

mm]- 



Before analyzing these equations let us indicate how to 
extract R(v) from its solution. Let us first rewrite the 
part of the action at the saddle point which depends ex- 
plicitly on V (in fact, on v ■ v): 



S[Xii,'^i, V]\expl 



m 

2T 

4 
m 



^i''Gf{k^O)v'' 



(248) 



= -^ni^G%k = 0)-GUk = 0)] 



[G^]{k,u) = Tgk- _^^ 

9k +[^i](u) 

Gi(fc,0)=g^rai(0), 
G25 = 9kTcr2i, 



and, to precision 0{n): 
Gl{k) 

G2i{k) 



Tgk~ (n/2)glTa2v, 
{n/2)glTa2ij. 



(240) 

(241) 
(242) 



(243) 



«2 



2T 



0'2v 2 

n 

2T 



Since the derivative of M4p[t'] = ^5'[x{j, tT{;,'D], w.r.t. v 
(in fact, w.r.t. v ■ v) only involves the explicit part we 
find, and comparing with (|230p : 



^|s'(4S2) = a2,. 



(249) 



Up to negligible corrections of order 0{n) the saddle- 
point equations for the diagonal blocks are the same as 
those in the regime v = 0(1), cf., (|14ll5p . independently 



We now discuss various cases according to whether or not 
the replica symmetry is broken in the diagonal block, (the 
saddle point of the GVM), indicating whether or not the 
system is in a glass phase. 
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1. Replica- symmetric region 



Assuming replica-symmetry, (j{u) = a , we have: 

G{k) = Tgkt + Ta^''9lh (250) 

G2,(fc)-Ta2^s<?^X (251) 

This reduces the saddle-point equation (|244p to: 

2 
cr^^ = --S'(2T/i), (252) 

independently of -D, and Eq. (|247p becomes: 

'x^^' = -|S' (4i^^ + 2T {/i + [.^s - ali]h]) ,(253) 



where /„ = Jf.g^- Comparing with (I252p we see that 
cr^l^o = '''RSj as is expected in the absence of external 
fields. 



2. Identification with FRG below the Larkin scale: 

It is convenient to introduce the function B' via the 
change of variables: 



^B\Av^)^a^!. 



(254) 



Plugging it into (|247p , we see that it satisfies exactly the 
self-consistency equation 



B\x) =B'(a; + 2r/i+4[B'(x)-i?'(0)]/2) (255) 



as derived in Ref. [l3| for the second cumulant function of 
B{x) which occurs in the effective action (defined there 
via R{u) = NB{v?)). Since (|255p uniquely specifies B 
up to a constant, a comparison of (I254p and (|249p shows 
that 



B'{x) = B'{x) 



(256) 



i.e., we have recovered, in the glassy regime, the general 
identity R = R announced in (|57l78p . 

From the self-consistency equation (|255p a FRG equa- 
tion for B{x) can be derived, following [cf., Eqs. (6.1)- 
(6.8)] in Ref. M 



mdmB'ix) = B"{x) 



'i{mdml2)[B'iO)^B'ix)] 



2mdmTIi 
l-h4J2B"(0) 



(257) 



with initial condition B{x) — B{x) for m = -|-oo. 

At this stage we only know that the FRG equation 
(I257p is valid in the replica symmetric region, i.e., for 
m > niciT) (see discussion in Section IIIBp . At the 



Larkin mass, 13" [0) diverges, signaling a cusp in the force 
correlator, as can be seen from the condition P?|) and 
(|255p . Hence at m = m+ the equation (|257p is still valid, 
but the second term (involving the temperature explic- 
itly) is negligibly small there and can be dropped. It is 
now crucial to establish how to continue this equation 
below the Larkin mass. To this end we compute B{x) ex- 
plicitly beyond the Larkin mass, i.e., in the glassy region, 
where RSB occurs within the replica groups, and derive 
the correct FRG equation for m < rric- 



3. Self-consistency equation for the second cumulant below 
the Larkin mass 



The saddle-point equation (|247p together with (|249p 
yields the correct continuation of the self-consistency 
equation (|255p below the Larkin mass: 



B'{x) = B' 



We expect that: 



2 / Gik)+AB'{x)l2 

Ik 



S'(0) = -|c7(u = 0), 



(258) 



(259) 



remains valid when replica symmetry is broken (with of 
course a different and non-trivial value for ct(0)) since 
it expresses the equality a{u = 0) = 172,5=0 between 
the self-energy associated with distant equilibrium states, 
and that associated with the coupling among the two 
replica groups in the limit of zero forcing, {i — > 0. In- 
deed, one can check using p4ip that ()259p is a solution 
of (pgg)) . given that ([Mi|) holds. Note that ((^55)) insures 
that the function B{rP') perfectly matches the large w^ 
limit of the v'^ — 0{1) results ()107|193p . irrespective of 
the scheme of RSB, as it should. 

Eq. (|258p can be rewritten in a form similar to (|255p : 



B'{x) = B' [x + 2Th + ^h[B'{x) - S'(0)] + 2 A 



where 



and 



B'(xo + x + Al2[B'{x)-B'{Q)] 



G(/c)-G(fc,0) 



Xo 



(260) 



(261) 



Ak - G{k) - Tgk - T<j{0)gl 



du[G(fc,u)-G(/c,0)] 



(262) 



is an anomaly which is non-zero if and only if the replica 
symmetry is broken. (To get the second line in (|262p . 
we used the generally valid inversion formula G{k, 0) ^ 
TaiO)gl) 
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4- Correlator and FRG equation below the Larkin mass 



From (|260p we easily obtain the behavior of the force 
correlator, —B'{x), at small argument. If the RSB so- 
lution is marginally stable with respect to a clustering 
instability at the largest scales, i.e., if the condition 



4/2S"(xo) 



(263) 



holds, with xo defined in (|261|) . one finds a non-analytic 
cusp in the correlator: 



B\x) = 



Ta{0) 



-2x 



(AhfB"'(xQ) 



0{x). (264) 



This type of marginality is automatically ensured in the 
case of continuous RSB, which is marginal with respect 
to fluctuations on all scales (see App. |l| including those 
corresponding to u^, which entails (|263[) . For a 1-step 
solution, however, it occurs naturally only on the tran- 
sition line Tc{m) for m > to*, while in the glass phase 
it requires the specific choice Uc = u'^p, cf. (|I85p . In the 
case of a non-marginal 1-step solution, the latter choice 
differs from the value of Uc corresponding to thermody- 
namic equilibrium, and was shown to describe rare disor- 
der configurations in Section UlID 21 On the other hand, 
for a fully stable saddle point with 1 > 4/2-B"(xo)7 one 
finds the regular correlator: 



-B'{x) = 



TaiO) 



xB"{xo) 
1 - 4/2B"(xo) 



+ 0{x^). (265) 



We emphasize the difference between these results and 
the generic non-analyticity (|147l205p found in the regime 
V = 0(1) at T = 0, occurring irrespectively of the RSB 
scheme. The cusp in the regime v = 0{\/N) analyzed 
here is present even at finite temperature, provided the 
RSB scheme is marginally stable in the sense of condition 
(I263p . The non-analyticity thus reflects the criticality of 
the system being at the brink of an instability towards 
an additional clustering of the topmost level of an ul- 
trametric structure. In the 1-step case, the instability is 
towards a two-step RSB with an additional step in the 
lower plateau22ii2^. 

Taking a derivative of (|258p with respect to x and using 
it to simphfy (|257p . one findsi^ repeating the steps of 
Eqs. (6.4)-(6.6) inRef.[ll 

- mdmB'ix) = B"{x) X (266) 

[imdmhlB'iO) - B'(a;)] + A(to)} , 

where 



A{m) = — 2m(9„ 
= —2md„ 



Th 



Ak 



4l2md,nB'iO) 



G{k) + 2Ta{0)mdml2 



(267) 



This equation is the general FRG equation valid in all 
regimes. In the non-glassy RS regime Aj- = and the 



amplitude A{m) is identical to the last term in (|257p . In 
the glassy region Eq. (|266p is certainly valid for all x > 0, 
and it is also valid at a; = if B"{0) < oo, i.e., if the RSB 
solution is not marginal in the sense of (|263p . But even if 
the solution is marginal, and hence 13" {0) is infinite (i.e., 
B'{x) exhibits a cusp), both sides of the equation have 
a non-trivial limit at a; = 0"*" which yields another valid 
equation as can be checked using (|264p . Another useful 
expression for A(m) is obtained by rewriting the second 
equation in formula (j267p as: 

A{m) = -mdraXo - 2Thmdm(T{Q) (268) 



where we used the definition (|26ip and (|262p . The second 
term in (|268p can be rewritten using the saddle-point 
equation (|244p at m = 0, and taking a derivative w.r.t. 
m one obtains 



A{m) 



-[l-4/2B"(xo)]ma„xo 
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a formula valid in all cases and regimes. We already point 
out and will discuss further below that the amplitude 
A{m) vanishes if the marginality condition (j263p is met. 

Let us now analyze in more detail the case of continu- 
ous RSB which includes as a limiting case the marginal 
1-step case occurring in d < 2 with 7 = 7c (d). In 
Ref. [13 it was found that the consistent FRG equationiSl 
for TO, < TOc(r) was (|266p without the last term, i.e., 
A{m) — 0. As we discussed above, this matches the 
established equation (|33p at to. = to+. The validity of 
(I266p for TO, < nic was inferred from the study of the 
FRG equation in inverted variables x = x{B'), noting 
that for models I and II the FRG flow completely stops 
at TO, = ?7i+ (the beta function is identically zero hence 
its natural continuation is to remain zero below rric)- In 
Section VIII-D of Ref. [H the analysis of the FRG flow 
was extended to an arbitrary bare model (i.e., an arbi- 
trary B{x)). It was found that the natural continuation 
is such that the function B' retains a cusp for all ni < rric- 
This was found to coincide with the marginality condi- 
tion at u,m, and hence lead to (|263l264p . 

From the above general expression (|269p we see that 
the vanishing of the amplitude A(m) = is a direct con- 
sequence of the marginality of the continuous RSB at 
u = u,,„, i.e., condition (|263p . and thus confirms the cor- 
rectness of the assumption made in Ref. 13. In addition 
we discover here that the cusp can be avoided if there is 
a non-zero amplitude A{m). This possibility was natu- 
rally not considered in Ref. 13, where the analysis was 
based on the self-consistency equation without anomaly 
and the ensuing FRG equation. While the amplitude 
A{m) always vanishes in the case of continuous RSB, a 
zero amplitude requires a specific choice of the break- 
point Uc = u'^P in the case of 1-step RSB, as discussed in 
the next Section. 
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Discussion of the FRG flow 



General considerations 



Let us analyze the generalized self-consistency 
Eq. (|260p for the renormalized correlator function B'{x). 
It depends only on the (given) bare disorder function, 
B'{x) and a number, xo = Xol'^j^^); defined in (|26ip . 
This is the only input from the saddle-point solution in 
the regime v ~ 0(1), which introduces an anomaly in 
the glass phase. We emphasize again that in the case of 
one-step RSB, xo does not only depend on T and tti, but 
also on the choice of the breakpoint Uc- 

Since the force correlator (the derivative of the poten- 
tial correlator) scales like B' ^ y^-2e+d+2C^ with typical 
values of the argument scaling as ra""^^ , it is natural to 
define the rescaled force correlator (as in Refs. Hlpa ): 
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At this stage, if we ignore all other information from the 
V ~ 0(1) regime, the two exponents C and 9 are unde- 
termined. Let us assume for now that they can be fixed 
by requiring that h' (x) reaches a non-trivial fixed point. 
From (|260p one deduces that h' satisfies the equation: 



b'ix) 
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where we have used I2 = Ad/ {tm'^) in the limit of infinite 
cut-off, and consequently set = d — 2 + 2^, as expected, 
so that the last two terms in the argument of B' (the 
bare disorder) scale the same way. The only information 
needed from the v ^ 0{1) regime is the value of xoj given 
by the saddle-point solution. It fixes the value of b' at 
the origin: 



6[,(m,r):=5'(0) = 
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(272) 
The above two equations then completely determine the 
FRG flow of 6', describing the evolution of the force cor- 
relator as the mass m is decreased. The corresponding 
flow equation was given in (|266p . Here we directly solve 
(|27ip focusing on models I and II, which admit simple so- 
lutions. We distinguish the case of continuous RSB (full 
or marginal one-step), which we only briefly recall since 
it was discussed in detail in Ref. Il3l . and non-marginal 
one-step RSB which requires a novel and thorough anal- 
ysis. The difference between the two cases can be grasped 
immediately. For continuous RSB, comparing (|272p and 
(PS)) . valid in that case for models I and II, we immedi- 
ately see that b'ifi) reaches a fixed-point value with the 
same choice of exponents as in the v ^ 0(1) regime. By 
contrast (|17ip and (|175p suggest a rapid decrease of 6'(0) 
as TO —^ as discussed below. 



Some features are independent of the RSB scheme and 
worth mentioning. Consider model I, B{x) — e~^ . If one 
makes the choice C = 0; then b{x) is the solution of 



b'{x) — SqCxp 



b'{x) - b'o 



(273) 



which has clearly a fixed-point form provided 5q, given 
by 



b' 



4.A, 



exp[-xo(™,T)], 



(274) 



reaches a fixed point as to, — > 0. Similarly for model II, 



B'ix) 



1 



, and (|272p together with the choice 



C — ((-y) ~ e/[2(l + 7)] leads to the equation: 
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Upon inversion this allows us to simplify the equation 
(|27ip for the rescaled force correlator as: 



m^K\^.(m^^ 



b'{x) - b'n 



(276) 

which reaches a fixed point if b^ does. Note that these 
expressions, as well as (I27ip . are valid for all to, both 
above and below the Larkin mass. 



2. FRG flow for continuous RSB 

As discussed in Ref. [ij for model I and II the flow is 
particularly simple. The structure of (|273p and (|276p im- 
ply that b'{x) = $(a; -|- xq) is always given by the same 
master function $, shifted by an amount xq which de- 
pends on TO only through 6q. For model I the master 
function is the solution of <I>(a;) = — eexp[— x— 1 — <I>(a;)/e]. 
The shift xq must be positive, and xo = leads to a cusp 
(infinite slope) at the origin. In the case of continuous 
RSB Xq freezes to zero as m decreases below mc(T). One 
checks that 6g reaches its fixed point value at m — mc{T): 



b'o = ^, 


(I) 


b', = e 


'^Ad 
e 



(277) 



(11) 



and remains constant everywhere in the glass phase. This 
implies in particular that 6' reaches its fixed point al- 
ready at the phase transition rndT) and sticks to it for 
all smaller m, a result inferred in Refs. KLlUlSl from the 
vanishing of the beta- function at mc{T). The FRG flow 
for models with continuous RSB is illustrated in Fig. [T] 
As shown in Ref. [ij for models different from I and II, 
and such that there is continuous RSB for all m < mc{T), 
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-b'(x) 




FIG. 7: FRG flow of the renormalized and rescaled corre- 
lator b{x) for model I in d = 3. The reduced temperature 
was chosen to be T = 0.5 and the curves correspond to 
1 < m/mc(T) < 5.5 in steps of 0.5, from top to bottom. 
For m < mc{T) the force correlator sticks to its fixed point 6', 
(thick Une) as given by (p73)) with 6'»(0) = e = l. Due to the 
specific form of B{x) in models I and II, the fiow reduces to a 
m-dependent horizontal shift of an otherwise constant master 
function $(a;). 
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FIG. 8: The shifts xo as a function of m at constant temper- 
ature r = 0.8 for model I in d = 1. The minimal shift xo =0 
is attained at the glass transition m = mc(T) — 0.673 and 
signals a cuspy correlator. In the glass phase, only the cusp 
condition u — u'^^ maintains :eo = and hence a cusp. Equi- 
librium and threshold states have regular correlators {xq > 0) 
in the glass phase. 



such as linear combinations of I and II with various pow- 
ers, the flow is slightly more involved. For m > mc{T) 
the flow is not a simple translation in x but the shape also 
changes as irrelevant parts oiB'{x) decay. At m = mc{T) 
a cusp appears and remains for all m < mc{T) as b'^ (and 
the shape of the function) converges slowly towards one 
of the above fixed-point values (|277p . 

In the continuous RSB case, the existence of this fam- 
ily of non-analytic fixed points is related to the vanish- 
ing of the amplitude A{m) in the FRG equation. In 
Ref. [ij they were derived directly from (|266p assuming 
A{m) = 0. For d > 2 and d < 2, 7 > 7c, they lead 
to the same exponents C, = ^(7) a-nd 6 = 0(7) as the 
study of the v"^ = 0(1) GVM regime, cf. Ref.iifi . Hence, 
all parts of the effective action scale consistently. Note 
that in the limit case d < 2, •y ~ •jc — 2/(2 — d), the 
glass phase is described by a marginal one-step RSB so- 
lution. As discussed in details in Ref. [l3| (Section VI. D 
and Appendix F), the high-temperature phase T > Tc is 
described by a line of analytic fixed points which solve 
(l?7T|) with C = (2 - d)/2 = C(7c), & = d{lc) = 0. In that 
phase b'{x) converges to one of these fixed points. This 
line of fixed points terminates on a cuspy fixed point, to- 
wards which b'{x) converge for T < Tc, and which has 
the form ((2761 [277)1 with 7 = 7c. 

By contrast, in the non-marginal one-step case, both 
for models I and II, the shift xq generically vanishes only 
at mc{T) and becomes positive again in the glass phase, 
unless the choice u = u'^p is made for the breakpoint, as 
we now discuss (see also Fig. |S]). 



D. FRG flow: new physics for one-step RSB 

1. Main analysis of FRG flow for one step RSB 

The situation is radically different in the case where the 
saddle point is solved by a non-marginal one-step RSB 
scheme. This happens for d < 2 and 7 > 7c (d) = 2/(2 — 
d). As discussed in Section IIIID 21 one may consider 
several possible choice for the breakpoint Uc. To interpret 
what is actually computed in each case, one should go 
back to the definition ([50)1 of the observable V[v] whose 
second cumulant is computed here in the regime v'^ ^ N. 

Let us first discuss the choice Uc = u^'^, where one 
computes the cumulant of the observable V[v] — V°'^[v] 
defined as the equilibrium free energy of a manifold (i.e., 
a directed polymer d = 1, or a particle d = 0), in an 
external quadratic well at position v and non-zero tem- 
perature T, in the limit of infinite N . More precisely, for 
every m this prescription selects the metastable states of 
the lowest free energy density in typical disorder. This 
condition is by definition equivalent to requiring the con- 
figurational entropy to remain fixed to S = as m varies, 
and thus closely resembles the "iso-complexity" condition 
often used to describe the rapid quench in systems un- 
dergoing one-step RSB^J^ 

The first observation is that for m < mc{T), A{m) > 
and this " anomalous" term acts as an " effective tempera- 
ture" in the FRG equation (|266|) which avoids the occur- 
rence of a cusp, similarly as a finite temperature does in 
the case of finite TV. Let us evaluate this amplitude at low 
temperature (i.e., away from the glass phase boundary. 
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T < Tc) and m ^ 0. One has: 
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and as shown in Appendix [Gj T/uJ'^ — l/u°'^ has a finite 
hmit as T — > 0. Hence xo diverges as ttt, — > and for 
model I B"{xq) = —B'{xo) = e^^° decays exponentiaUy 
fast to zero. For model II, AhB" {xo) ~ TO-2+(2-d)7^ 
hence in all cases where 7 > 7c one finds: 
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This term in the FRG equation (|266p has precisely the 
same form and scaling as a standard (one loop) temper- 
ature terroiS^ with T -^ T/uc, avoiding the occurrence 
of a cusp. Remarkably, this temperature takes in the 
limit 771 ^ the same value as the effective temperature 
Toff which arises in the modified fluctuation dissipation 
relation in such systems^^. 

Let us now discuss in more detail the solution of the 
self-consistent equation (I27ip . We start with model II. If 
in (|270p we choose the exponent values from the region 
v^ - 0(1), i.e., C = (2 - d)/2 and 6 = d- 2 + 2C = we 
find that: 
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where < for 7 > -fc- Hence, with this scaling the 
second cumulant in the region v"^ ^ 0{N) tends to zero. 
To study this flow for tti — > one can instead use the 
choice 9 ^ 6 in (|270p and find a fixed point for the scaled 
function defined in that way. Indeed, one can check easily 
that with the same C = (2 — d)/2, as to ^ 0: 
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which indicates that the force fluctuations are reduced 
by a factor m~^^ in the region v^ ~ 0{N) as compared 
to the region v^ ^ 0{\). For model I, also obtained by 
taking 7 -^ 00, the effect is even stronger, the decay of 
correlations being exponential of the form: 



B'{x) 
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which cannot be put in the form of a scaling function as 
in (|270p . The equilibrium flow for that case in d = 1 
is shown in Fig. [9l It simply amounts to the horizon- 
tal shifting of the master function $(a;) by Xq^, which is 
plotted in Fig. [51 While the correlator first grows with 
decreasing mass in the high temperature phase (like in 
d > 2), it decreases again within the glassy phase. The 
reason is that as m decreases below nic, the stability of 
the one-step energy landscape increases and suppresses 
sample-to-sample fiuctuations of the global shift u° of the 
states, which in turn sets the scale for force-fiuctuations. 




FIG. 9: FRG flow of the equilibrium correlator b{x) for model 
I in d = 1, at the reduced temperature is T — Tnmx/2. 
The dashed lines (light blue) correspond to m/mc(T) = 
{2, 1.5, 1.25}, from bottom to top, and the dotted lines (pur- 
ple) to m/m^{T) = {0.4,0.3,0.2} from top to bottom. Note 
that the correlator only exhibits a cusp exactly at the contin- 
uous transition, m = mc{T), with b'Q{mc{T)) = e — 3 (thick 
line). Due to the specific properties of model I the correla- 
tor always equals the same master function, displaced by a 
m-dependent horizontal shift xq. 



2. Force correlator, shocks and the cusp 

In the limit 7 ^ 7^ one has 61 — > 0, as one recovers 
the marginal one step solutioufii^. and B'{0) then scales 
in the same way in the two regimes v'^ ^ 1 and v'^ ^ N. 
As noted previously the value of B' (0) from the v'^ ^ N 
regime described by FRG is proportional to the contri- 
bution of the lower plateau to the two point correlation 
at fc = 0. It also gives the sample to sample fluctuations 
of the global shift of states, u^, in a given sample p63p . 
While in the case of continuous RSB they scale as the 
relative fluctuations of the states, this is not the case for 
the non-marginal one step, since there Go{k) [and uo] is 
much smaller than G'i(fc) — Go{k). In the non-marginal 
one step case it is the upper plateau which dominates the 
two point correlation (for any k ^ 1/nic). Upon moving 
the potential well by w '--' N'^^'^, the global shift of states 



becomes a function of v, u^ — u1{v), the connected cor- 
relator of M°(0) and u^{v) being described by G2v{k)- 
In this regime the free energies are of order TV, and one 
may expect that configurations centered around u^ (0) or 
u1{v) differ typically by quantities of order E{N) which 
diverge as A^ ^ 00 (at least as N^^"^). In that case, when 
two levels cross as a function of v, the switching of the 
equilibrium position takes place as a function of the scal- 
ing variable E{N)/T. Hence, it turns into sharp jumps 
as A — > 00 even at finite T . 

If these jumps are of order N^/'^ in the displacement 
w, they show up as a non-analyticity in the force corre- 
lator. This is what is observed in the case of marginal 
one step and continuous RSB. On the other hand, stable 
(i.e., non marginal) one step systems seem to have a much 
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smoother landscape leading to an analytic response of the 
system to a displacement of the well. This indicates that 
if there are shocks, their size does not scale with N-^^^, 
i.e., the system proceeds in small, and effectively rounded 
jumps in contrast to the marginal cases. 

The non-analytic response of systems with continuous 
RSB is tightly bound to their marginal stability with re- 
spect to clustering. This criticality of the system im- 
plies an anomalous response to a field (v) pulling apart 
two replica groups, whereby the above argument suggests 
that this proceeds via the occurrence of shocks of size 
N^^^. Notice that the presence of the non-analyticity 
only relies on the marginality at Um whereas it is insen- 
sitive to the marginality of the RSB scheme at larger u 
(deeper down in the ultrametric hierarchy). This indi- 
cates that the shocks are really associated with jumps 
on the largest scales involving u"(w). That is, the whole 
hierarchy of states is affected, and not only part of it, as 
it can happen in the regime ii^ ~ 1. 



3. Some subtleties of the FRG flow for one step RSB 

Although we have focused on the flow at small tti, in- 
teresting features also happen around the Larkin mass. 
At 771c a cusp appears and A{m) first vanishes at m = mc, 
but then becomes positive again. This behavior indicates 
that a bifurcation in the lowest free energy state occurs 
at rric, where metastable states appear and the system is 
critical. However, the further evolution of the system for 
m < rric is smooth in the sense that the system becomes 
uncritical again, with no further bifurcations occurring. 
This is seen to translate into a response which is regular 
when the harmonic well is displaced by v with v'^ ^ N. 

To complete the discussion of the equilibrium FRG 
flow, we mention a subtlety (corresponding to the choice 
Uc = Uc**) which arises for temperatures in the range 
T G [Tc{0),T^''^]. An inspection of the phase diagram 
in Fig. [4] shows that as m — > 0, the thermodynamic equi- 
librium is given by a replica symmetric saddle point in 
the GVM, which is reached from the glassy phase by a 
discontinuous transition upon decreasing T. In that tem- 
perature regime, one natural continuation of the FRG 
flow from large to small mass is to maintain the con- 
dition d(j)/duc = 0, defining u^'"^" even though it cor- 
responds to u^'^'^ > 1 in the region where the genuine 
equilibrium is described by a RS solution. This reflects 
the fact that the selected states are thermodynamically 
irrelevant, in the sense that their Gibbs weight is neg- 
ligible. Nevertheless, there are small corners of phase 
space where metastable states exist even when the mass 
is small, and among those states, the ones selected by 
u^'^'i" still correspond to those with the lowest available 
free energy density in typical disorder. We have verified 
that this branch of one-step solutions indeed exists down 
to m = for ah temperatures T < T""^"^. 

We now turn to the discussion of other possible choices 
for the breakpoint. If one chooses Uc = u**^ one presum- 



ably computes the cumulant of an observable, V[v] — 
V[vY^, defined as the free energy of the threshold states 
for the manifold in an external quadratic well at position 
V. These threshold states being known to be relevant for 
the dynamics, we may conjecture that the quantity com- 
puted using this criterion will find an interpretation as 
a part of the dynamical effective action. Note that un- 
der this prescription the configurational entropy of the 
selected states, I]*'^(m) grows with decreasing mass, re- 
flecting either the bifurcation or birth of threshold states. 

Another choice of breakpoint, and the only one which 
leads to a cusp, is Uc — u'^p. Then the amplitude 
A{m) — 0. Hence it leads to similar fixed points as for 
the continuous RSB case. In particular, it yields a fixed 
point of 6g and freezes the renormalized correlator to the 
cuspy shape it attains at mc(T). The exponents are then 
the same, i.e., C — C(7)- Note then that 9 = 9{'y) < 
which means that the temperature is relevant. That such 
a choice can be made was conjectured in Refs. IlllllSl for 
T — 0. Here we show it to be possible at any tempera- 
ture. It is not difficult to prove that one-step solutions 
satisfying the cusp condition (|185p exist for all T < Tmax 
and TO < mc{T). Note, however, that one finds u'^p > 1 
for m < m^ and T > Tc{m), where Tc{m) is the (un- 
physical) branch of the instability line defined by (|177p . 
see Fig. [5] 

A natural question is then: what observable would be 
selected by this process? As discussed in Sec. 11111X2) the 
choice u'^P appears to select metastable states that only 
exist in rare disorder configurations (the smaller m the 
rarer). It remains to be understood whether such an 
observable could be constructed by imposing some con- 
straint on the metastable states. 



E. Non-uniform v 

The above analysis is easily generalized to a non- 
uniform Vx, and allows to compute the functional R[v\ for 
a non-uniform Vx in the regime v^ ^ N where it takes 
the form R[v\ = NB[v'^] where ^[v^] := B[{vl}x], i.e. 
it is a functional of the field v"^. We just give the final 
saddle-point equations and the resulting self-consistent 
equation for _B, which generalize the one obtained in 
RefjiS (for R[v\) in the replica symmetric region (see 
Eq. (3.31) there with the correspondence in notations 

s(v'^(xyvt^(x)) ~^ <^v{x)). One easily sees that the diago- 
nal blocks (in the space of two replica groups) are again 
independent of Vx and therefore independent of x, and 
thus simply correspond to the MP solution. Generaliz- 
ing (|247p . the saddle-point equation in the off-diagonal 
sector yields a single equation, which now involves space 
indices: 
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where 
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which generalizes (|242l) . As in the uniform case 
o'2.{i=o(2;) = o'(O). Solving for the function a2,v{x) yields 
the desired functional since: 
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which is the analog of Eq. (|249p . These equations gen- 
eralize (j260|) to a non- uniform field v^- It would be in- 
teresting to see whether FRG equations can be derived 
also for the non-local parts of the functional B, but this 
is left for future investigation. 



V. DISCUSSION AND CONCLUSION 

We have computed the renormalized disorder- 
correlator R{v) in the large N limit for an elastic mani- 
fold in a 0(N) symmetric random potential in presence of 
an external quadratic well. It was obtained as a physical 
observable describing the correlation in free energy when 
the center of the quadratic well is varied by v. It contains 
direct information about shocks, i.e., abrupt switches be- 
tween two competing equilibrium positions that occur as 
V is varied. We have demonstrated the existence of two 
large- A^ scaling regimes v'^ ^ 1 and v^ ^ N and obtain 
closed expressions in each regime, and their zero and low 
temperature limits. Our results provide a direct connec- 
tion between the GVM approach [v = 0) using replica 
symmetry breaking saddle points and previous large A^ 
FRG approaches. 

We found that as v is increased shocks start to oc- 
cur at a uniform displacement of the well of order v ^ 
V* ~ L~ " which decreases with system size or, equiva- 
lently for a non uniform displacement v^ of order one but 
confined to a bounded region of volume one in space^^^. 
These shocks are rounded by temperature, but turn into 
a non-analytic cusp of the force correlator at T = 0. Our 
results bear some similarities to shocks found in Burgers 
turbulence in large dimension. The shocks can be inter- 
preted within the RSB picture of an ultrametric phase 
space, which predicts that in a given disorder environ- 
ment there exist several states, all centered within a sin- 
gle big valley whose position itself fluctuates from sam- 
ple to sample. The shock regime w ^ w* then describes 
the energy crossing between states as the harmonic well 
is shifted, and accordingly, it is sensitive to the replica 
symmetry breaking structure of the GVM approach. In 
particular, one finds that all GVM saddle points related 
by replica permutations contribute to the computation 
of observables. Beyond this regime, the force correla- 
tor remains constant in a large interval N^^^ 3> w > f * 
which reflects the sample-to sample fluctuations of the 
force density associated with the big valley, itself related 



to the sample-to-sample fluctuations of the global dis- 
placement w° of the valley. 

This value of the force correlator matches perfectly 
with the value obtained from FRG calculations in the 
regime v'^ ~ A^. In that regime it was found that FRG 
recovers only the fluctuations among most distant states, 
which is governed by the (non-trivial) lower plateau of the 
self energy function, (Tq ■ It is now clear that the FRG in 
this regime captures the shock structure of the big valley 
position u'^{v). Indeed, imposing a position of the well on 
the scale v^ ^ N affects the global shift vP itself, which 
may lead to shocks. We have computed the decaying 
correlations between m°(0) and u^{v), as well as the force 
correlator. In the case of systems with continuous or 
marginal one step RSB, such as manifolds in 2 < d < 4 
or directed polymers in long range correlated disorder, 
we have found that their marginality towards a cluster- 
ing instability implies a cusp also in this second scaling 
regime w^ ~ A^. This reflects the fact that the shift of the 
center of the well provokes shocks which are abrupt even 
at finite T because they involve the crossing of energies 
which scale with A^. By establishing the persistence of 
the nonanalytic force correlator throughout the regime 
of small m and T, we have proved that the previously 
obtained FRG flow had been correctly extended into the 
glassy regime. 

A major step forward has been achieved in the case 
where the regime v^ ^ \ displays a non-marginal one- 
step RSB, a problem which was left open in the previ- 
ous FRG study. In particular we have found that in the 
regime v^ ^ N the force correlator is non-analytic only 
at the glass transition, but not within the glassy regime. 
We interpret the latter in terms of a smooth energy land- 
scape where the manifold evolves smoothly (on the scale 
N^/"^) as the center of the harmonic well is varied. In the 
glassy regime the disorder correlator becomes strongly 
suppressed with decreasing mass which is a consequence 
of the increasing stability of the one-step landscape and 
the associated smallness of sample-to-sample fluctuations 
of vP . Further, we have shown how the exact FRG equa- 
tions in the regime w^ ~ A can be extended into the 
glassy regime, and how replica symmetry breaking trans- 
lates into an anomaly in the FRG flow equation. More- 
over, the possibility of studying non-equilibrium branches 
of nietastable states in a one-step system by tuning the 
break point parameter Uc was found to be equivalent to 
tuning the anomaly in the FRG equation. 

The present study raises many questions. First, at 
finite A^ the two scaling regimes in v are not clearly dis- 
tinct. The main task is then to determine which aspect 
of either regime will persist at finite A^. Until now the 
FRG flow equation in the loop expansion seemed to con- 
nect more directly to the regime v^ ^ N . In particular 
one then expects that the shocks of the regime v'^ ^ N 
should become rounded by temperature at finite A^, as 
was found in FRG calculations in the loop expansion and 
for A = 1. However, the present study raises the possibil- 
ity of a new regime in the FRG accessible only at smaller 
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V. The interpretation of the cusp at T = in terms 
of shocks remains vahd, and offers an interesting venue 
for future studies going beyond the mean field picture of 
switching between states. These concepts also generalize 
to other disordered systems, e.g., to spin glasses, where 
the statistics and properties of shocks can be studied by 
similar methods. Work is in progress in this direction. 
This should yield valuable new insight into the structure 
of the phase space and elementary excitations, such as 
droplets. 

Note added in proof: After completion of this paper we 
became aware of the work of H. Yoshino and T. Rizzo^, 
which studies related spin models and obtains results on 
shocks and cusps similar to those discussed in the present 
work for the regime ii^ ~ 1. 
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Transverse fluctuations: 



G(fc = 0) 
G(fc), X 



0G(fc = O), 



1 



'/ 



G(fc), X 



Self-energy and other auxiliary functions: 



[cr. Si,. g ^{k/fic 




R{v) 

R[v],Rn,'Rn 



In, 
f \ 'f 

2n 



r 
1 



^^iE^l 



R{v/rf), 



E3 



R[v/rf],Rn,'Rn 



W,q,Q,u,h -^ W,q,Q,\i,b. 



(Al) 



With these substitutions all equations turn into dimen- 
sionless identities, whereby in some cases additional fac- 
tors of c = ii^~^Ecr^r need to be restored. Note that 
amplitudes such as A being dimensionlee are unchanged. 



APPENDIX B: DIRECT EXPANSION OF W TO 
ORDER V* 



APPENDIX A: RESTORING DIMENSIONS 

In the main body of the text we used the natural units 
r/ for transverse lengths, fj,c for the mass (and l//ic for 
longitudinal lengths in d > 0), and Ec for energies. In 
order to recover the full dependence on the parameters 
-Bq, rf and c, one simply has to restore the dimensionful 
units so as to render all quantities dimcnsionless: 

Masses and longitudinal lengthscales: 



TO, fc, A — !■ — [to, fc. A] 



Temperature and energy densities: 



T -> T/E,, 



The perturbative expansion of (j6ip to second order 
requires 



r,(2) 
abed 



'-'ab '^cd ■ 



(Bl) 



In the case of two sets of replica with v°- = v^'^/2 for 
a = 1, ..n/2 and u° = —v^'^/2 for a = n/2 + 1, ..n, one 
has Gab — lab/'^ with qab defined in ([57)1 . As for the sum 

-fih' ■= L^ <^ib^ = '^^ab + P (cf., EHl), rephca symmetry 
restricts this tensor to take the form 

E'abcd — ^^abcd + B{Sabc + Sabd + Sacd + ^bcd) 

+ CiSabScd + C2{5ac5bd + Sad^bc) + Di{Sab + Scd) 
+D2{Sac + Sad + Sbc + Sbd)+E, (B2) 

the coefficients of which can be obtained by solving the 
system of linear constraints arising from the identities 

V^ p(2) _ ^cp(l) 
/ . ^abcd — '-' ^cd ' 



5 bP^^^ 



bed 



GP, 



(1) 

d ' 



^ 5ac5bdPabcd ~ X! ^^"-b 
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The quantity of interest for the second order expansion 
of (ICTI) is the second cumulant: 



w^H(2) = 



abed \abcd / 



(B3) 

For the case of two rephca groups, one obtains the final 
result: 

VK[t;](2) = ^4 [„^ + „2(^^ ^ 2C2) - n^a^] (B4) 



V rr{2 — n) 
(3-n)(l-n)2 






2^4 

~3~ 



1 / 1-1 

2/ 



G^(u)du - / G(u)rfu 
\Jo 



7l2 + 0(„3)^ 



From this one derives (|114p in the text, by recalling that 



42^2 
i?2H = lim ^^iyM(2) 



n^O 71^ 



(B5) 



and G(u) = g(u)/2. 



APPENDIX C: A USEFUL IDENTITY 

For any function <I>(j/) with derivatives decreasing 
strictly faster^^^ than l/\y\ for large \y\, one has 



dyy[^{y + z) + $(y - z) - 2^{y)] = 

-z2[$(oo)-$(-oo)]. (CI) 



Indeed, 



/oo 
dyy[^{y + z) + $(y - z) - 2^{y)] 
-00 

dyj; / dz'[$'(y + z')-*'(y-^')] 



dy / dz'[^{y + z') - $(y - z')] 

-00 JO 

dy I dz' I dz"^'{y + z") 
= -[$(cx))-$(-oo)] / dz' [ dz" 

Jq J-z' 

== -z2[$(oo) - $(-00)]. 

Applied to integrals over Gaussian averages, normalized 
s.t. (z^) = 1, Eq. ini]) implies 



(C2) 



dyy^{'S>{y+y/Qz)),-<^{y) 

-|[<i>M-$(-CX3)]. 



APPENDIX D: PERTURBATION EXPANSION 

FOR L'^v'^ < 1 



Here we compute the two lowest orders in the expan- 
sion defined in Section llll C 21 The functions mi^2 satisfy, 
from (fTTTjl : 



mi = --g(u) [m'^ + u{mly] , 

"12 = -^q{^) ["T-i + u(2r7iomi)'] , 
To lowest order, using that: 

< + {mlY = 0, 
one has to solve: 

TOi =--(l-u)g(u)m^'. 



(Dl) 

(D2) 

(D3) 
(D4) 



The solution is: 

1 f" 
mi=p(u)m[; , p{u) = - du(l -□)<?(□). (D5) 

Integration by parts gives 2p(u,„) — J^ du (1 — u)q{u) = 

— g(0)-|-/p q{u)du. The integration range can be extended 
to the interval [0, 1] since q — outside [Um, Uc]. Plugging 
this into (|112p and using J^ yrriQ = — J^ TOq = — 1, the 
term proportional to q{0) cancels and one obtains (|113[) 
in the text. 

The next-order correction satisfies, from (jDSP and 
(ID2]): 



7712 = --(?(u)p(u) [77lo" + U (27710777,0 )' 



(D6) 



This gives: 
1 



'712 = ■;^"io" / du qp + {momQ)' j duuqp. (D7) 



In order to calculate (|112p . we need 

/oo 
dyym'^" = -[77io]!!°oo = 0, 
-00 



(D8) 



dy77io(r7i^)' = -, (D9) 



dy y{m^'m'^)' 



where on the second line we used (jD3p . Inserting into 
(fTT^ gives 

R2[v] = ^T^ j duuqiu)p{u), (DIO) 

where j7(u) is defined in (jD5p . Substituting the definition 
(|D4ID5p . ug = 2p+q, integrating by parts, using p(l) = 
one finds: 

/ du uq{u)p{u) = [p2 + qp]l - - / du (u - l)qq 
Jo ^ Jo 

= lj\uq^- \q{Of - p(0)[p(0) + g(0)], (DU) 

which with p(0) = 1/2 [/„ duq — ^(0)] yields the final 
result (I114P in the text. 
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APPENDIX E: THERMAL BOUNDARY LAYER 

OF THE EFFECTIVE POTENTIAL 

CORRELATOR 

Here we revisit and complete the calculation of the ef- 
fective potential defined in Ref. [43 , in an attempt to 
connect with the FRG function R{u). The effective po- 
tential studied there is constructed (see Ref. l43 for de- 
tails) from the probability distribution for a given Fourier 
mode u — Uk (denoted there 0o) in a given environment. 
This sample dependent probability is denoted here Zy [u] 
(unnormalized) and there Vn{4>o)- The authors then in- 
troduce the potential correlator: 



■n/2- 



V{u ~ u) — — lim 



Z(u)"/2^(u')"/2 _ z{u) Z{u') 



rn/2 



{n/2YP 



2/92 



(El) 

describing the second moment of the effective potential 
V{u) = (3~^\n[Z{u)]. This is the central object stud- 
ied there. Although it is a very physical object, it is 
very different from the effective potential studied here 
V{v) - which is also a physical observable - but involves 
a source (implemented via a quadratic well). As a re- 
sult V{u) cannot be connected to the standard FRG^^'^. 
Some discussion of the differences between the two ap- 
proaches was given in Section VIII G of Ref. 13. Here 
we discuss more of the differences in details, based on 
explicit calculation. 

Since this version of an effective disorder is also inter- 
esting, it is worth to push here further the calculation 
of Ref. |43 . To avoid confusion between u (displace- 
ment field) and u (replica overlap) we replace everywhere 
M — > w, keeping in mind however that its physical mean- 
ing is different from the one (position of well center) given 
in the text. The above partition function can be evalu- 
ated via a saddle point as in the GVM: 



n^(t''^)-E^^p -f^^^^^G-^^w^"^'^ 



o=l 



exp 



''^(fc2-|-m2)(w2+t,'2) 



4 
^exp 



2 



[v-^^^aatv-C'^) 



(E2) 



In the same way as we derived R[v], we can compute 
the correlator of V as: 



V{v - v') = -2T2 / dy y M(0, y) - t&nh{y) 



where M satisfies the same flow equation as M, but with 
the " inverse coupling" : 



9(u) = -(w-w')^cr(u). 



(E3) 



We can now apply our general expressions (|126ll34p to 
the above case and find for the thermal boundary layer: 

- V{v - v') = Vo{v ~v') + Vi{v~v') + ... , (E4) 



where: 

Vo{v - v') 
Vi(u-w') 



2T^g(u,) 



Ta{uc) 



{v~v')\ 



(E5) 



T' 



dq 
du 



ijiz^2iq,-q)) 



= (u-i;') UcCrc(l+7) X 

^ J f 

dvvo ( ^p \ z 

-^ \ v 



V - 



T 






where Gc = cr(uc). These formulae have well-defined 
limits as T ^ (recalling that u^/T = Am^^ and 
Ta,^ = 2m2-7[(2 - e)A]). For m -^ the last term 
reduces to; 



Vi{v 



{v - v'fu,a,{Ta,Y/^ r[l + g/(2-g)] 
4 r[5/2-t-6'/(2-0)] 

With the notation introduced in Ref. |43 (for d>2) 

A-d 



1+7 



?-l 



rf-2' 



1 



9 ■= 



-1 



(3a{uc)ul{v - V 



l\2 



(E6) 
(E7) 



the above results can be recast into the form: 



V{v- 



where: 



= -Vo - Vi = 



(/3uc 



1+7 



- U T2 ^ 



T2(..) ^ ^ 



(..... JMi^ 



(5^8) 



and 4'{x) = 2a;coth(x), as before. Using ■0(x)^2|x| for 
large arguments, one further finds: 



T2(X) 



y2a;3/2 



1 



r(i/(i + 7)) 



(1 + 7)5/2 r(5/2 + 1/(1 + 7) 



The term in brackets is easily seen to tend to 1 as 1 + 7 = 
e/(d - 2) ^ 0. Thus, T%-^°{x ^ 00) ^ y^^^/^^ 

For large arguments w, the flow of M is attracted to 
an intermediate fixed point M(u,y) « tanh(uy). This 
happens for a large coupling parameter g ^ 1 (or more 
explicitly for displacements such that {v — w')^77i^+^ ^ 
1). If in addition the mass is sufficiently small, [y — 
v'Y ^ m^'^+^\ the correlator V is controlled by this 
intermediate fixed point and can be shown to scale as 

V{v - v') ^ g2e/(2+e) _ (^ _ ^>Ye/(2+0)_ (^g^ 

An important difference to the FRG correlator R{v) 
concerns the scale on which the nonanalyticity lives in 
V{v). Comparing Vo and Vi we find that higher order 
terms become dominant for v > w** ^ rric , reflect- 

ing the fact that the effective potential correlator V is 
sensitive to physics on the Larkin scale, rather than to 
shocks occurring on the scale to, as is the case for the 
FRG correlator. 
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APPENDIX F: COMPARISON WITH THE TBL 
FOR DROPLETS 

Due to the 0{N) symmetry, the disorder correlator is 
only a function of v = |u| (for spatially uniform v). Ac- 
cordingly, the force correlator splits into transverse and 
longitudinal parts: 



-d.uA,Riy) = ^L{v) 



V 



iv){ 



% - ^ ) ' (Fl) 



R'iv) 



with Al{v) — —R"{v) and At{v) — — " ''^•' , derivatives 
being w.r.t v = \v\. Both correlators are equal at i; = 0, 
so that dv^dv TZ(v) is well defined and proportional to Sij 
at V = 0. Note that at large N one has j^^id^JZ = 
At(w) + 0(l/iV), and hence the transverse components 
dominate the average force correlator. 

The TBL contribution to the droplet force correlator 
()155|) can be expressed in terms of longitudinal and trans- 
verse correlators: 



Ad/°P(.)=T^(y?V 



— ^ — V 



-A-P (^) - T^ ivl^ C^v 



(F2) 
(F3) 



where j/i = j/ii denotes the component of y parallel to v, 
while j/2 is an arbitrary orthogonal component. 

To make contact with expressions obtained for the 
manifolds, we recast these expressions into the form: 

- Aijiv) =T J dbpt:r''ib)mv\). (F4) 

with the distributions 
Pr^ib) = 4 fd^-'y^b'D fyii = ^,yA , (F5) 



m^ J 



N-i I m'^y2\ „/ 26 



In computing the force correlators (|159p we have used 
that: 

d{, (v ilriiixz)) = ( 2v'(l){vxz) -\- vz — ip(vxz) ) (F6) 
\ dz I ^ 

= (i)(l + z^)->p{vxz)) ^ , 

5? (v^ipivxz))^ = (F7) 

(1 + z^^-tplvxz) + z{l + z'^)-—ip{vxz) 

dz 

{{z^-z'Wvxz))^. 



for any d and A. Equations (|166I167I168I171I172I175() 
form a closed set of equations, which together with 
(|170|179p determines Uc ^ uf^: 



= ^^ = ^[Bixo] B{xi) (xo - Xi)B'ixo)] 



T 



9k 



ln(l + Sigfc) 



.(Gl) 



The easiest way to show this is to consider Si and Uc as 
independent variational parameters of the saddle-point 
equations, s.t. dSi/dUc — 0, and then use p76p . 

Let us now analyze the limit T ^ 0+ of the vari- 
ous one-step saddle points for rf < 2 and A = oo. It 
is easy to see that the corresponding saddle-point equa- 
tions (I175I176|) together with the three possible condi- 
tions ^7^, PM)1 or P^ admit solutions of the form 
Uc = Tuc with finite u^ as T — > 0: 



4.A, 



1 



Xo 



1 



1 



-,2-d 



e(2 - d) Uc 

Ei = -2Q,[S'(xi)-S'(xo) 



Si)i-'i/2 



(G2) 



-2uc[B'{0)-B'{xo 



The three possible values for the breakpoint are: 

(i) existence of a cusp in the FRG in the large-ii regime 

(Uc = u^P), i.e. condition (fT55)) : 



^i?"(X0) - 1, 
em'^ 



(G3) 



(ii) the equilibrium static condition (uc = u°'^), which 
from (jGip using ()173p reads: 



A, ( 



B{0)-xoB'{xo)] 



(G4) 



Si 



e{2-d)ui 



(m2+Si)i- 



(m2+Si)2 -m" 



and (iii) the condition (|184p for threshold states (uc 
u'^), which is equivalent to: 



El = ml{T) 



mliO) 



(G5) 



For concreteness, we give explicit results for model I in 
the limit to — > 0; an illustration for general m and d = 1 
at a fixed temperature can be found in Fig. [Sj 



For small m, B{xo) ^ B'{xo) ~ e 



Xo 



and (To -^ 0. 



One finds, from dGTj) for any < T < Tc: 
T 11^=^ 



^c ^ — rp : ^1 



^^^-dulVi2-d) (gg) 



APPENDIX G: 1-STEP SADDLE POINTS IN THE 

GVM AND THEIR LOW TEMPERATURE 

LIMITS 

Here we first give for completeness the general one-step 
saddle-point equations for the equilibrium statics, valid 



where Tc is given in (|187p . Using condition (iii), one finds 
for any < T < T^: 



thg2(l-u-)/. ^ T 



A=<{T) 



Td 
th 

2— p" 
T^ 



-2ul^/i2^d) 



(G7) 
(G8) 
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where Td is the dynamical transition temperature given 
in HM]). In the r ^ hmit one obtains>ii^ 



de 



d=l 



Tcexp{d/e) 



0.7937, 



-.th 



2 2 I e 



2/e 



d=ll 

2' 



(G9) 



The cusp condition (|G3p imposes that xo ~ logm, and 
from the saddle-point equations it follows that: 



u^P-m''-71og(l/m) 



(GIO) 



We finally determine the T — > limit of the config- 
urational entropy associated to this branch of cuspy 
solutions, S(u^P) = u^^. For m <C nic one finds 
Si « 2u^P ^ m. Hence, the leading term in the con- 
figurational entropy is: 



E(u^P) « (Q?P) 



2 d4> T— O.m^O (U^P) 



cp^2 



dilr 



Adin 



d-2 



d)elog(4^d/eTO') 



(Gil) 



which is negative. Hence, the condition (|G3[) selects ex- 
ponentially rare configurations. 



APPENDIX H: THERMAL BOUNDARY LAYER 
IN THE CASE OF ONE-STEP RSB 

Here we prove that the thermal boundary layer in the 
1-step case can be written as 



7^lH =2T^u,v(^ 



where 



VWQ) 



(HI) 



dQ'{z^j2Q'coth.{z^j2Q')),,. (H2) 



Since both expressions at small Q are easily seen to be 
of order 0{Q), it is sufficient to show the equivalence for 
the second derivative with respect to y/Q. Taking the 
second derivative of TZi[v] in (|200p . we find 
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^7^l 



(H3) 



2T^u,d{VQ)^ 

J —oo 

-2^$" ($^' - $) - Zz'$^$i, - 1 



where $ = $(?/) = lncosh(?/), $^ = $(y + zy/Q), and 
primes denote derivatives with respect to y. After a par- 
tial integration (in y) of the terms containing $", the 



expression simplifies to 



dy(z(z-z') (!-*>'.')).,.' 
'(z-z')2 



-V'((z-z')v^)) 

/ z,z' 



z^ijiz^)) 



(H4) 



where V'(o) = 2acoth(a) as in ()133p . and we have used 
the fact that z + z' and z — z' are independent Gaussian 
variables with variance 2. 
On the other hand, one finds 



1/'(VQ) = 2v/0(zV2gcoth(zv/2Q))., 
V"i^/Q) = (4zv/2Qcoth(zv/2Q)). 

+ (2v/2g22a^ coth(zv/2Q))^, 

= (2y2Qz^coth(zv/2g))2 

= (zV(2v/2Q))., (H5) 



where we performed a partial integration in the second 
but last line (with the remaining term coming form the 
derivative of the measure e~^ ''^), completing the proof. 



APPENDIX I: GENERAL SOLUTION FOR 
CONTINUOUS RSB 

The general continuous RSB solution has been derived 
in detail for the model g^^ — k^ + rn^ in Ref.[l3. Here we 
sketch its generalization to arbitrary g^ . A continuous 
RSB Ansatz is always just marginally stable on all scales, 
as expressed in the present case by the identity 



1 = 4B" 2 / G{k)-G{k,u] 



, (11) 



. [.9fei + H(u)]^ 



for all Um < u < Uc. Using ^T5\i . which in continuous form 
reads a{u) = -|S' (2 J^ G{k) ~ G{k, u)\ , this leads to 
the relation 

a(u) = --B' {{B")-' { ^ -X\ ■ 

^ V V4/,k7^+M(u)]-v; 

(12) 
Further, the relation 1/u = da/d[a\ allows one to solve 
for cr(u). 

In particular, using [(7](um) = 0, we find immediately 
the general expression for a{um) ~ f (0): 



aiO)=a{u^)^-^B'({B")-'(^^ 



(13) 
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as the limit of (G(ui) 

-G(u,_i)]l^,_) + [G(l)-G(uK-i); 

G(uif)) of a K-step RSB as A" — > oo. 
u_ft-<..<ui<n being reversed when the 
analytical continuation to n < 1 is performed. 
The replicon instability and the clustering instability are 
defined and discussed in detail in Ref. 75, cf., Eqs. (A3. 3, 
A3.7-8) and Eqs. (5.11, A3. 16-17), respectively. 
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Other choices of interest are, e.g., Oi = ——!f \-r-, — to for 

Of < 2 which models long-range elasticity. 

The bare correlator B is defined in the same way as in 

Refs. [lllll3l . but corresponds to —f = B defined in Ref. l46l. 

Corresponding to B{z) = \b{0) - ^1 9(1 - 7) -h 

^Bo (1 + ^) ^ for 7/1, and B{z) = B(0) - 

Bo \n{l -\- z/rj) for 7 = 1. 

In Refs. 1 1111131 the model I was parametrized with g = Bo, 
rf = 1. The bare correlator of model II was taken 
in the form B'{z) = — g , ^j)^^„ , corresponding to the 
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reparametrization ga '^'^ — Bo/r^, a^ — jr^- 
For finite A'" there is a range of values 7 > 7* (A*') for 
which model II renormalizes to the SR disorder class I, 
with ^*{N) —> 00 as N ^ 00. 

One should distinguish this "Larkin scale", which is a 
simple combination of dimensionful quantities meant only 
to provide an order of magnitude estimate, and a "Larkin 
mass" or a "Larkin length" which in some cases can be 
defined precisely from an observable, e.g., for A*' = 00 an 
exact phase transition to a regime with many pure states 
occurs at a mass m = mc, see below. 

In the case of RSB (i.e., in the glass phase and within the 
GVM) there are some subtleties in this correspondence. 
One can define the intravalley connected average, given 
by Gaa — G(l) (see below) - relevant for time scales in 
which no transitions between valleys occur - as well as 
the full connected correlation Gaa — /g duG(u). 
In the limit of large UV cutoff A 3> mc{T), corresponding 
to the weak-collective pinning regime. Note that the UV 
cutoff needs to be retained only in the T dependence of 
m^{T) via f, cf., (fT9l). 

According to (|19|) . in dimensions d > 2 the temperature 
scales as T ~ f\-(.d-2) ^j^j^ f < q^^-^s^ ^^ ^^^ transi- 
tion and in the glassy regime. Thus, in the weak collec- 
tive pinning regime where A 3> /ic = 1, one always has 
Uc ~ T ~ A"'''"^' < 1. In long-range models d < 2 with 
7 < 7c, one finds that on the transition line the fluc- 
tuation dissipation ratio Uc = (1 -I- 7) (2 — d)/(4 — d) < 
(1 -I- 7c(d))(2 — d)/(4 — d) = 1 is always smaller than 1. 
The definition of V[v] can alternatively be interpreted as 
the free energy of the elastic manifold (with elastic prop- 
erties described by gk), where the quenched disorder has 
been displaced by —v{x). 

As usual formulae for discrete sums are obtained replacing 
J^ = Y2x^ Ik = ^ '^ X]fe s-iid for a uniform v^ ~ v one has 
Vk = L'^vSkfi- 

The term —qi-i naturally arises in the construction^^ for 
a complete Parisi matrix, see also (|65p . Here it arises from 
the extra term —qe-i in the definition of gi which in- 
volves an incomplete Parisi matrix, (the replica sum is 
interrupted at u^). 

When comparing (|115|l with (11131) . note that G^^o = 
because of HTl) . 

Note that the function T^ri[{)] in (|134p diff'ers from 
T'^fi[v] in the systematic expansion (|115|l by a trivial 
quadratic term -2T^dq{uc)/dT\T=o oc T^v'^ = Tv^. 
Independently for each of the A*' components u''^ . This is 
left implicit here in order not to burden the notations. 
Equivalently, one can write Wa = exp{—f3Fa) with a tail 
P{F) ~ exp(uc/3-F) for large negative F. 
We make the natural assumption that B"{x) is 
monotonously decreasing with x. 

We correct a misprint of the 1-step solution given in 
Ref.m , which suggested that u(r ^ 0) ~ T^2-d)/{4-d) ^ 
However, u**" will exceed 1 at small m and T close enough 
to Tc{m,). 

Even though distinct they may occur at the same value of 
parameters in the case of the marginal one step solution. 
In the case of lattice glasses one finds the same regions of 
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stability for Uc as here: The 1-step solution is unstable to 
replicon fluctuations for Uc < u"^, while it is unstable to 
the clustering instability for Uc > u'^^. 
As in Sec. IIIIC3I it does not exactly coincide with the 
expansion in powers of T since some extra (weak) T de- 
pendence is hidden in other parameters such as, e.g., Ei. 
However, as before the series can easily be rearranged 
later. 

Agreement with (|147p and (|149p can be checked for the 
marginal one step case using 6* = 0, Uc = ^ from (|23p. 
and m^ -I- Ei = m^. 

Note that (m" - u^f = [fifi]"" - 2[fi{i]"'' + [uu]''''. 
Note that the condition is on v and not on v since the 
saddle-point equations are expressed in terms of C. In 
other words, the condition is i;""'"' = v"^ + o{\/7<f). 
One may also consider higher bare cumulants S^''' . While 
both the saddle-point equation in the regime u ~ 1 and 
the self-consistent equation for the correlator are changed, 
the essential results about the connection between the 
FRG regime v^ ^^ N and the RSB regime u ~ 1 remain 
the same: One can derive a closed self-consistency equa- 
tion for the renormalized force correlator B' (involving the 
bare cumulants S''''), from which the exact FRG equa- 
tion can be derived; in the glass phase RSB introduces 
an anomaly in the selfconsistency and the FRG equation; 
the identity B'{0) = -Ta{0)/2 matching FRG and RSB 
regime holds everywhere in the phase diagram; the force 
correlator exhibits a cusp at the origin, if and only if the 
saddle-point solution is marginally stable with respect to 
clustering fluctuations. 

For B, but this should be the same as for B. 
In 1-step systems, the complexity usually jumps to a flnite 
value at the transition, which is then assumed to remain 
constant throughout the glass phase under conditions of 
a rapid quench, because the trapping threshold states be- 
come stable and do not bifurcate. In the present case the 
situation is similar, except that the complexity vanishes 
at the transition. 

^°^ To one loop a non-zero temperature generally produces an 
additional term ~ TVlR{u) in the FRG flow for R{u). 
Upon proper rescaling of R{u), so that a, T = Q fixed 
point exists, its coefficient generally scales as ~ Tm and 
temperature is said to be irrelevant if ^ < 0, relevant if 
6^ > or marginal if 6^ = 0. 

^^^ The X dependence in the limit is delicate, since (|282p is 
only correct asymptotically as r7i — ^ if 7 > 7c, but not 
for 7 = 7c. This is because, for 7 = 7c, one cannot neglect 
the last term in the argument of B' in (|271|) . while it is 
negligible for 7 > 7c as m — » 0. 

^^^ such as a suitable norm, e.g. J v'^, remains of order one 

^^^ If the function is even, l/\y\ behavior is allowed. 

^^"^ Although it is not excluded that another version of FRG 
could be built from it - in a spirit close to Polchinsky RG 
see discussion in Ref. Il3l . 

^^^ This result is smaller than the one given in Ref. |43 by a 
factor of \/2tv. 

"5 In Ref. m the result Q;;'' = l/rc(m = 0) is given. This 
results, however, from a rather unrealistic assumption on 
the behavior of B(x) at small x. 



